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Metastability in a condensing zero-range process in the 

thermodynamic limit 

Ines ArmendMz* Stefan Grosskinskyt Michail Loulakis^ 


Abstract 

Zero-range processes with decreasing jump rates are known to exhibit condensation, where a finite 
fraction of all particles concentrates on a single lattice site when the total density exceeds a critical 
value. We study such a process on a one-dimensional lattice with periodic boundary conditions 
in the thermodynamic limit with fixed, super-critical particle density. We show that the process 
exhibits metastability with respect to the condensate location, i.e. the suitably accelerated process 
of the rescaled location converges to a limiting Markov process on the unit torus. This process has 
stationary, independent increments and the rates are characterized hy the scaling limit of capacities 
of a single random walker on the lattice. Our result extends previous work for fixed lattices and 
diverging density in [J. Beltran, C. Landim, Probab. Theory Related Fields, 152(3-4):781-807, 2012], 
and we follow the martingale approach developed there and in subsequent publications. Besides 
additional technical difficulties in estimating error bounds for transition rates, the thermodynamic 
limit requires new estimates for equilibration towards a suitably defined distribution in metastable 
wells, corresponding to a typical set of configurations with a particular condensate location. The 
total exit rates from individual wells turn out to diverge in the limit, which requires an intermediate 
regularization step using the symmetries of the process and the regularity of the limit generator. 
Another important novel contrihution is a coupling construction to provide a uniform bound on the 
exit rates from metastable wells, which is of a general nature and can be adapted to other models. 
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1 Introduction 


A rigorous understanding of metastability phenomena in stochastic particle or spin systems has been 
a subject of major recent research. Intuitively, in such systems, the configuration space contains two 
or more disjoint metastable sets (called wells in the following) with an associated separation of time 
scales phenomenon in the scaling limit of a system parameter. In each well, the process spends a 
very long time which allows it to equilibrate to a metastable state. Exits from wells are then triggered 
by rare fluctuations, which lead to exponentially distributed waiting times in a well. Once activated, 
transitions to other wells occur on a much faster time scale, and do not depend on the detailed past 
of the sample path. So the limiting metastable dynamics corresponds to an effective Markov process 
on a highly reduced state space of metastable states associated to the wells. On a mathematically 
rigorous level, important conceptual questions of current research interest include a proper definition of 
metastable states in terms of probability measures, as well as a general practical framework to establish 
the separation of time scales. An important aspect in this context is the precise definition of metastable 
wells and in particular an optimal choice of their size, which is to some extent arbitrary. Intuitively, 
maximizing the depth and minimizing the complexity, or size, of the wells leads to a stronger separation 
of time scales. In combination they effectively characterize the free energy landscape for the chosen 
definition of wells in analogy to the classical framework of competition between energy (depth) and 
entropy (complexity). 

Different approaches to metastability are summarized in ll40l . Chapter 4, including a pathwise 
treatment 1201 which is based on an analysis of empirical averages along typical trajectories, and has 
more recently been studied also in l28l . For reversible systems a powerful potential theoretic approach 
has been developed uniiii, using systematically the concept of capacities to establish sharp estimates 
on expected transition times between metastable states. This technique was applied in various models 
and is summarized in the new monograph ifTSl and the review papers |[l4]|22l. Potential theoretic 
methods have been applied to a particular family of condensing zero-range processes in lll|6l, leading 
to the development of a martingale approach summarized in 19] which we follow in this paper. Instead 
of deriving exponential limit distributions of individual exit rates from metastable wells directly, the 
limit process is identified as a Markov process via the solution to a martingale problem. 

The zero-range process was introduced in ll44l as a stochastic particle system without restriction 
on the local occupation numbers, where the jump rates of particles depend only on the occupation 
number of the departure site. The dynamics locally conserves the number of particles and the process 
is known to have a family of stationary product measures m, which can be indexed by the particle 
density p. Under certain conditions on the jump rates this family has a maximal element at a finite 
(critical) density pc, and the process exhibits a condensation transition. Conditioned on densities p > Pc 
the system phase separates into a homogeneous phase distributed at the critical product measure, and 
a condensate, where the remaining mass concentrates in a single lattice site for typical stationary 
configurations of large, finite systems. This phenomenon was first reported in the theoretical physics 
literature Il25ll27]l . and rigorous results including the equivalence of ensembles, a law of large numbers 
and a central limit theorem for the condensate size have been established in ll^lT5l 1^. In spatially 
inhomogeneous zero-range processes a condensation phenomenon can occur that has different dynamic 
and static features as discussed above, and can be studied with the help of coupling techniques (see 
e.g. ||4| and references therein). In this paper we focus on spatially homogeneous systems, which are 
necessarily non-monotone BD and where basic coupling techniques cannot be applied. 

On a translation invariant lattice of size L the condensed state exhibits an L-fold degeneracy with 
respect to the location of the condensate, which is uniformly distributed under the stationary measure. 
The metastable limit dynamics of the condensate location for reversible processes has been established 
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in GlISl in a slightly more general setting allowing for spatial inhomogeneities, but these results are 
restricted to a fixed lattice in the limit of diverging particle number N. In this setting, the depth of 
metastable wells dominates their complexity and they are relatively small sets. Repeated visits to 
particular configurations occur before the process exits the well (often called attractor states), which 
provides a renewal structure that can be used to establish the separation of time scales. The limit process 
of the condensate location is then a continuous-time random walk on the fixed lattice, where transition 
rates are proportional to the capacities of a single random walker. For the same model, matching upper 
and lower bounds for exit rates from metastable wells have been derived in m in the limit L,N —i' oo 
with diverging particle density N/L —> oo (see also ifTSll . Chapter 21). 

In this paper we establish the full metastable limit dynamics for the condensate location for the 
above family of reversible zero-range processes in the thermodynamic limit, i.e. we take L,N ^ oo 
with finite, supercritical density N/L p > pc- We consider a one-dimensional system with 
periodic boundary conditions, which leads to a Levy-type limit process for the rescaled condensate 
location on the unit torus with stationary, independent increments. A major new difficulty is related 
to thermalization within the metastable wells, which are exponentially large in the system size and 
previously applied renewal techniques do not apply. To separate the slow and fast variables effectively, 
we describe the thermalization dynamics on the wells using the restricted process introduced in IT2l 
and Q. We can then prove relaxation time bounds for this process by a comparison argument with 
independent birth-death chains. Another key ingredient is a uniform bound on exit rates from a well 
which is established by a novel coupling construction and is used in several occasions in our proof, 
including thermalization. The characterization of the limiting generator of the process is more difficult 
in the thermodynamic limit due to increasing complexity of the free energy landscape, caused by the 
diverging number of wells and transition paths between them. To tackle this we have to introduce 
an intermediate regularization step on a coarse-grained lattice using the full symmetries of the model 
and regularity of the limit process. We note that, similarly to the results in ifT^ for diverging particle 
density, we are not able to establish matching bounds for transition rates between individual wells. But 
in order to identify the limit process it is sufficient to get matching bounds only between regularized 
sets of wells. As in previous results on this model, the leading order of transition rates between wells 
is polynomial in the system size L, which requires sharp estimates on transition rates from potential 
theory. Due to the increasing complexity of transition paths we need a better control than results in 
|I7]|3 on the leading order error terms, which constitutes an additional technical difficulty. We note that 
a rescaling argument in a different context has also been used in ll38l to establish metastability for the 
ABC model. 

So far, there are only few metastability results that deal with infinite volume limits as summarized 
in Ga, Part VII. Examples include kinetic Ising models at low temperatures M or low magnetic 
fields ll43l . results on the dynamics of critical droplets ll3?l . or dilute gases at vanishing density 
and temperature ll^ . All these results feature the additional scaling of a system parameter such as 
temperature or density, which increases the depth of metastable wells compared to their complexity. 
Often this is a necessary requirement for metastability to occur at all, and in lfT9l this has been used 
as a simplifying assumption to obtain results for the zero-range process in an infinite density limit. To 
our knowledge, the thermodynamic limit result derived here without any scaling of system parameters 
constitutes therefore one of the first metastability results of this kind. 

The paper is organized as follows. In Section|^we introduce notation and state our main result. The 
proof and a more technical discussion of the result is given in Section The remaining sections are 
devoted to establishing the main ingredients of the proof, including the equilibration dynamics in wells 
in Section]^ uniform bounds on exit rates in Section]^ and the derivation of the generator of the limiting 
process in Section]^ Upper and lower bounds on expected transition rates from capacity estimates are 
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given in Section]^ and Sections]^ and [^contain auxiliary results on tightness, martingale convergence 
and properties of the stationary measure. 

2 Notation and main result 

2.1 Notation 

Consider the zero-range process {ri{t) : f > O) on the one-dimensional discrete torus A = Z/LZ with 
N particles and state space 

Xl,n = {p G Nq : S'l(? 7 ) = N} where we denote S'l(p) = '^Vx ■ (2.1) 

xGA 

The dynamics is defined by the generator 

^ P{x,y)g{r]:„){f{r]^'y)- fir])) ( 2 . 2 ) 

A 

for all continuous functions / : tv —> R, with the usual notation = rjz — ^z,x + ^z,y, ^ G A 

for the configuration 77 ^’^ where one particle has moved from site x to site y. We focus on symmetric, 
nearest neighbour jump probabilities p{x, y) = -f ^Sy^^+i with periodic boundary conditions, 

and jump rates are of the form 

( 77 - \ ^ 

-) for n >2 , with parameter 5 > 0 . (2.3) 

n — 1/ 

Without the canonical constraint Sl = N, the process is known to exhibit a family of stationary product 
measures with a maximal element 1 / |II]l27l[36l|35l, which has marginals 

11 " 

iy[r]x = n] = -—— with gUn) = TT g(k) = for n > 1 and gUO) = 1 . (2.4) 

ZcgKn) 

As long as 6 > 2 the normalization and first moment of 1 / are both finite, i.e. 

CXD ^ CXD 

Zc := I + ^ n~^ < 00 and pc := v{r]x) = — ^ < 00 . (2.5) 

n —1 n—1 

The process is irreducible on X/, tv and the corresponding unique stationary measures p,L,N are called 
canonical distributions, and can be written as conditional product measures 

gL,N = y'[ ■ I'S'i = X] . ( 2 . 6 ) 

To simplify notation we will write p = pl,n from Section [^onwards. 

We study the large scale behaviour of the process in the thermodynamic limit with particle density 
p > 0 , i.e. 

L,N = Nl —>■ 00 , such that N/L —>■ p . (2.7) 

For b > 2 and p > pc the process is known to exhibit a condensation transition in the following 

sense. Denoting by 

Mi = max Pa; (2.8) 

x^A 
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the maximum occupation number as a relevant order parameter, we have a weak law of large numbers 

Ell [361 Ea 


Ml! L —>■ 


0 , P< Pc 

P ~ Pc J P ^ Pc 


as L,iV —>■ oo (2.7| , 


(2.9) 


where convergence holds in distribution w.r.t. the sequence p.L,N- A corresponding CLT has been es¬ 
tablished in 0, and the fluctuations are Gaussian and of order \/L for b > 3, and of order i) 

with an associated stable law for 2 < 6 < 3. Moreover, the distribution of the configuration outside the 


maximum is known to converge to the maximal product measure v with density pc in the limit (2.7 1 , 
so the largest occupation number in the bulk outside the maximum is typically of order 
which holds for all 6 > 2 independently of the fluctuations of the maximum. Therefore for supercritical 
densities p > p^ configurations exhibit a unique extensive maximum with high probability, called the 
condensate, and the rescaled location of the maximum 


’ipLi'n) = ^{x e A : r]x = ML{r])} (2.10) 

is given by a single site in the rescaled lattice ^A. By translation invariance, tp^is distributed uniformly 
in under p,L,N, and ergodicity of the process on implies that iphipit)) visits the whole lattice 
on a long time scale. It is expected that this dynamics is metastable, i.e. tfjL{r]{t)) is constant for a long 
(random) time interval, and then changes abruptly to a new value which depends only on the current 
location of the condensate. 


Our main result is that, indeed, for large enough b > 2 and p > pc the zero-range process ( 77 (f) : 
f > 0 ) in the thermodynamic limit (2.71 exhibits metastability with respect to the observable tpL on the 
time scale 


er. := L 


l+b 


( 2 . 11 ) 


This means that the sequence of processes {ptBc ) • ^ ^ O) converges to a Markovian limit process 

on the unit torus T = K/Z, the scaling limit of ^A. A rigorous version of this result is provided in 


the next subsection in Theorem 2.2 including an exact formulation of the mode of convergence and 
required assumptions. 


2.2 The trace process and main result 

The rigorous formulation and proof of our main result on the large scale metastable dynamics of the 
condensate follows the martingale approach developed in lEliaii, which requires a partition of the state 
space as a first step. We define the well as the set of configurations where the condensate is located 
at X S A, 


■— {77 e Xl,n ■ Vx = Ml > N - pcL - ul, % < I3l for all y ^ x} . 
We choose the scales in this definition as 

and /3 l=2[L5^J, 


aL = L2 + 2i, 


( 2 . 12 ) 


(2.13) 


which allow for typical fluctuations of the condensate size of order ^/L, and of bulk oc cupa tion numbers 


of order for any b > 3. Note that under our conditions on b in Theorem |2.2| we also have 

^ L. We choose 13l as a sequence of integers for later convenience, while the exponent in a/, 

(2.14) 


is optimal for the estimates in Section [6(21 We denote by 

E := VJx&kE' 


and A := Xl,n \ E 
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the set of all wells and its complement. The first result states that the process started from a well spends 
only a negligible amount of time on the set A on the timescale 0^. As usual, we denote by and 
the path measure and expectation of the process ri{-) with generator (2.2 1 and with initial condition 
77 ( 0 ) = p. 


Proposition 2.1 Replacement by the trace process 

For the process defined in \2.2\ with rates (2.51 and b > 20, we have for all t > 0 

rBr.t 


supE^ 

r]eS 


j 1a(0(s)) (is —>■ 0 as L, N ^ 00 {2.7) with p > Pc . 


(2.15) 


The proof of this result is given in Section 5.3 The larger lower bound on the parameter 6 is a 
purely technical restriction which we discuss in Section 3.2 Since the process spends only negligible 
time outside S, re-parametrizing time by the local time on S will not change the process in the limit. 
Denote by 


Tt'■= f l£{r]s)ds and iSj := sup {s > 0 : 7^ < f} 

Jo 


(2.16) 


the local time on £ and its generalized inverse, respectively. The trace process is then defined as 

: f > 0 ) with := pst , (2.17) 

which takes values in the set £ and is well defined since 7* diverges P^ — a.s. as t ^ 00 due to 
the irreducibility on the finite state space As is shown in Q, Section 6.1, (2.17 1 is in fact an 

irreducible Markov process on £ with jump rates 


■= r{'n,C) + X! '^ivX)^dTe = 7?] • (2.18) 

C6A 

Here r{ri,^) := are the jump rates of the process rj{-) with 

generator \2.2) , and 

:= inf {f > 0 : 77t G J-} (2.19) 

denotes the hitting time of a set J" C Xl^m- For point subsets 7^ = {^} we use the shorthand = T^. 
The trace process has generator 

= [fid - fid] , ( 2 - 20 ) 

and reversible invariant measure pl,n [ • |^] restricted to £. To simplify notation, we will call 

pd-]=FL,N[- \£]- ( 2 . 21 ) 

On the set of wells £ the rescaled location of the maximum can be written as 

V’l(p) ^ ^ a; l£-(?7) G T = K/Z , (2.22) 

x^A 


and is well defined without degeneracy. We can now state our main result. 
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Theorem 2.2 Consider the zero-range process with generator \2.2\ and rates (0) with b > 20. Fix 
a seq uence of initial conditions 77 ( 0 ) G £ such that —>■ Iq G T. In the thermodynamic limit 

{2.7) with N/L p > pc we have for the trace process {2.17) on the time scale 9 l = 




(2.23) 


where convergence holds weakly on the path space £*([0, cnd), T) in the usual Skorohod topology. 

The process {Yt : t > 0) has stationary and independent increments on T, initial condition Yq and 
generator 

= [ r'^{v) {f{u + v)- f{u)) dv 

for all Lipschitz continuous functions / : T —>■ K, with rates 


\v) = 


V gT . 


(2.24) 


(2.25) 


Zch {p - PcY^^ driO^v) 

Embedding the unit torus in [0,1) C K, the distance on T is given by dT(0, v) = [tiKl — |t7|). Zc is the 
normalization of the invariant measure 0, and 

is a constant depending only on b. 


(2.26) 


Note that (2.15 1 and our main result do not apply for initial conditions 77 ( 0 ) G A, which are. 


however, untypical if the process starts from the stationary measure pl,n (cf Corollary |9.5| l. The 
mode of convergence in terms of the trace process as presented here has been introduced in ||21|i^, and 
extended recently in 19] to a more general context. The (random) time change in the definition of the 
trace process, which is negligible in the limit L -G 00 , can also be absorbed in a definition of a suitable 
topology on the path space of the limiting process. Further discussion, including possible extensions of 


our result, is provided in Section 3.2 


3 Proof of the main result 


The proof of Theorem 2.2 uses a standard approach to establish existence of a limit process by a 
tightness argument, and to identify the limit by the solution of a martingale problem. We follow the 
method outlined in Q, proofs of auxiliary results used below are given in the rest of the paper. We also 
discuss the main novelties and possible extensions. 


Here and in the following sections we adopt a few shorthands and conventions to avoid an overload 
of notation. We write p = pl,n for the invariant measure of the full process, and for the invariant 
measure of the trace process ( |2.21[ ). Constants denoted by C are independent of L and N, and can vary 
from line to line. 


3.1 Proof of Theorem 12.21 

The proof of convergence holds on arbitrarily large compact time intervals for the limit process, and 
throughout this section we denote the length of this interval by T > 0. Let 

(y,^:f>0)e79([0,T],T) with Y/^ := ^L{vH0Lt)) , (3.1) 


8 












be the speeded up process of the rescaled maximum location, and let be its distribution on path 
space I?([0,T],T). 


Proposition 3.1 Tightness _ 

Under the conditions of Theorem 2.2 with & > 20 the sequence of path space distributions is tight 

onD{[0,T],T). 

The proof is given in Sectionj^where a control of the quadratic variation excludes the accumulation 
of jumps and ensures that limit points have right-continuous paths. Tightness implies existence of sub- 
sequential weak limits of in the Skorohod topology, and we denote any such weak limit by Q. In 
order to identify the limit we need to show that for all f < T and all Lipschitz-continuous functions 
/ e Lip(T) 


/(wt) - /(wo) - [ C^fiujs)d. 
Jo 


s is a martingale 


(3.2) 


where ojt '■ f?([0, c»), T) —> T is the coordinate process on path space. Toge ther with the uniqueness 
result for the martingale problem associated with proved in Subsection |8.3[ this implies convergence 
of and characterizes the limit Q as the law of the Markov process (l) : f G [0, T]) with generator 
(2.24 1 , because Lipschitz functions form a core for this generator. Since T > 0 is arbitrary, this 


implies Theorem 2.2 Precisely, we need to show that 

g((uJu : 0 < M < s)) (/(wt) - f(uJs) - C^f{uJu) du^ 


= 0 , 


for all 0 < s < f < T and all bounded, continuous functions g : Z?([0,T],T) 
t G [0, T]) is a Markov process with generator , we know that 


(3.3) 

1. Since {r]^{9Lt) : 


f{Y,^) - fiY,^) - 9 l f C^if o i^L){9^{9Ls)) ds = 

JO 

= f{Y,^) - f{Y,^) - C^f{Y,^) ds + - OLC^if o ^i)(/(0is))) ds 


is a martingale for all t G [0, T] and L G N. We will establish below that 


(£^f{Y,^)-9LC^{fo^^){g^{9Ls))) ds 
as L —i' 00 , which implies that 


supE^ 

rieE 




g{{oJu '■ 0 < u < s 


- f{ujs) - / C^fiuJu) du 


0 


(3.4) 


(3.5) 


for all 0 < s < f < T and bounded, continuous g : I?([0,T],T) —^ K. To identify the limit of the 
left-hand side with (3.3 i we use the following result, which is immediate from Lemma 8.2 in Section 


Proposition 3.2 Let mI(uj) = /(wt) — f(uJs) — f* du. Then 

S'® g((uJu ■ 0 < u < s)) m/(uj) —>■ E® g{{uJu : 0 < m < s)) M/ (uj) 
as L —)• oo for any t G [0, T]. 


(3.6) 
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To conclude the proof we have to show (3.4i, i.e. that we can replace the generator of the trace 


process with that of the limit process (cf Section 3 in Q). This is the main part of the paper and 
is divided into several steps. Since we cannot compare the generators and directly, we will 
introduce an auxiliary processes on T with generator which is explained in detail below, and we 


rewrite the time integral in (3.4 1 as 


o ds = 

= - 9L£^f{Y,^)) ds 

+9lJ^ ->C^(/oV'l)(/(0ls))) ds . 


(3.7) 

(3.8) 


The goal is now to show that the terms in (3.7 1 and (3.8 1 vanish individually as L —oo in L^-norm. 


In (|3.8|l we compare the trace process to the auxiliary process on the rescaled lattice C T, which 


describes the effective jumps of the condensate location. Its generator is dehned as 

C^f{x/L) := r^{z){f{{x + z)/L) - f{x/L)) , 

zGA,z/0 

with jump rates (using translation invariance) 


(3.9) 


(3.10) 






given by the expected rate between wells for the trace process. Therefore, we can also write the generator 
(3.9 1 as the expectation £^/(x/L) = p,|£x . Before the location of the condensate changes, 

the process remains in the same well for a long enough time to equilibrate, and the transition between 


wells is effectively described by stationary averages of jump rates as in (3.101, which is established in 
the next result. 


Proposition 3.3 Equilibration in the wells 

Under the conditions of Theorem \2.2\ with b > 20, 

rt 


SUpE^ 

rie£ 


9 l / (/:^/(n^)-£^(/o^i)(p^(0is)))ds 


(3.11) 


as L ^ oo, for all t G [0, T], 


In addition to using Lemma [33] given below, the proof requires an estimate on the relaxation and 
mixing times within a well, which have to be strictly smaller than 0/,, and is provided in Section]^ 

In ( |3.7| i we replace the generator of the auxiliary process with that of the limit process using the 
following result. 

Proposition 3.4 Dynamics between wells 


Under the conditions of Theorem 2.2 with b > 5. 


supE^ 

rjeS 


[c'^f{Y,^)-9LC^f{Y,^))ds 


-G 0 


(3.12) 


as L ^ oo, for all t G [0, T]. 
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The proof is given in Section and requires sharp bounds on the transition rates of the auxiliary 
process, which are provided by capacity estimates in Section In order to get matching upper and 
lower bounds, an important new step in this part of the proof is to regularize the rates r^{z) on an 
intermediate scale and use the regularity of the test f uncti on /, which is explained in detail in Section 
1^ This finishes the proof of our main result. Theorem 2.2 □ 


An important estimate that is used in the proof of tightness, equilibration and replacement of the 
trace process is the following uniform bound on the exit rate from a well. The proof of this Lemma is 
given by a coupling argument in Section]^ which is one of the crucial new results of this paper. 

Lemma 3.5 Uniform bound on the exit rate 

Under the conditions of Theorem \2.2\ with b > 20, there exists a constant C > 0 such that the exit rate 
from a well is uniformly bounded by 


sup V < 


C 


r]eS° 




L® log L 


(3.13) 


3.2 Discussion 

Main new ideas. Our proof follows the martingale approach outlined in ii, which was previously 
applied to zero-range processes on a fixed lattice in the limit N —>■ oo Q. In contrast to this case, the 


thermodynamic limit (2.7 1 considered here involves a significant change in the complexity of metastable 


wells, since the sizes of the wells and the number of transition paths between them increase with L. 
Following the discussion in 1291 , this presents a technically more challenging metastability scenario, 
in particular since the free energy barriers of the metastable wells in the zero-range process are only 
logarithmic. We quickly summarize the resulting conceptual and technical difficulties and the three 
main novel contributions of this paper to overcome them. 

• In previous work with limits of diverging particle density mini the depth of the metastable wells 
was dominating their size, and they could effectively be replaced by individual configurations, 
so-called attractor states. The repeated visits to those configurations before hitting another well 
lead to a relatively simple renewal-type proof for equilibration in the wells. In our case of finite 
particle density the size of metastable wells increases exponentially with L and these types of 
arguments do not apply. Instead, we have to describe the metastable states as probability distribu¬ 
tions on wells. As outlined in |l9l|29l a suitable candidate is simply the restriction of the stationary 
measure, and to prove Proposition |3.3| we use suitable dynamics restricted to a well following Q, 
see Section]^ Using a standard jacknife estimate ll26l we establish a Poincare inequality compar¬ 


ing these dynamics to independent birth-death chains in Sec tion 4.2 and obtain a bound on the 
relaxation time on the metastable well of order in Lemma 4.1 This is clearly not optimal, and 
a sharp bound is expected to be of order L^, which would moderately improve our conditions on 


the parameter b in Propositionand Lemma 3.5 currently 6 > 20, to & > 13. 

Resulting from capacity bounds presented in Section]^ the transition rates between wells r^(z) 


(3.10i scale like l/(|z|(L — jzl)) (6.6i. The total exit rate from a well on the scale 9]^ diverges 


as log L ( |6.7| i in the thermodynamic limit, and the exit time distribution from a well degenerates 
in the scaling limit and does not converge to an exponential random variable. This is due to small 
jumps which accumulate at diverging rate, but are still negligible on the macroscopic scale for 
the rescaled limit dynamics, which is established rigorously by showing tightness in Section 
Related to this, the errors in ( 6 . 61 also do not allow to get matching upper and lower bounds on 
transition rates between individual wells, in analogy with results in ||T91 for diverging density. We 
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address both issues in the standard thermodynamic limit for fixed density, using the symmetry of 
the system and the regularity of the limit dynamics on the unit torus, which are fully characterized 
by Lipschitz-continuous test functions. We regularize the dynamics on an intermediate scale in 
Section|^to get matching bounds on rates between regularized sets of wells, which are sufficient 
to derive the limiting generator and prove Proposition |3.4| As an additional technical difficulty 
we have to keep track of corrections of capacity bounds to leading order in L, which require 
very precise estimates on the stationary measure summarized in Section]^ Sections]^ an d0 are 
independent of the rest of the proof, and the only restrictions on the parameter b arising there are 
given by equations (|6.4|i and (|6.13|l, resulting in a much weaker condition of 6 > 5. 


While the li miting metastable dynamics are determined by stationary averages of transition rates 
given in (3.101, a uniform control of the exit rates from a well is important to estimate error 
terms as is discussed in general in ||9l. Some particular examples of spin systems where this has 
been achieved are mentioned in ifTSll . but note that using capacities between individual configura¬ 
tions in the thermodynamic limit would lead to bounds that diverge exponentially in the sy stem 


size. Here we derive a uniform upper bound on exit rates scaling as L ^ log ^ Lin Lemma|3.5 


which is proved in Section]^ by a novel coupling construction with a growing number of birth- 
death chains. The number of chains increases only linearly in time, and the construction ensures 
that in the event of changing well, at least one of the chains has grown a condensate. This is a cen¬ 
tral auxiliary result of the paper, and is used in the proof of equilibration in the wells (Prop.[33]l, 


replacement by the trace process (Prop. 2.1 1 , as well as for tightness (Prop. 3.11. The proof of 
Lemma [T5] requires b > 20, and sharper estimates on equilibration times in wells would slightly 
improve this condition as mentioned above. 


The bottleneck of the method leading to restrictive conditions on the parameter b results from the 
inherent competition between depth and complexity of metastable wells. Both quantities increase with 
the size of the wells, and the aim to maximize depth and minimize complexity leads to an optimal 
choice of their size that enters most prominently in the uniform bounds on exit rates. The coupling 
argument in Section [5T| gives a lower bound on the expected time to change wells for the full zero-range 
process. To turn this into an estimate for the trace process, we have to bound the time spent outside 
the wells on the set A, which is controlled by the invariant measure in Corollary |9.5| Larger wells of 
increased depth improve this bound, but at the same time lead to an increase in the number of transition 


paths to other wells. Both effects compete and affect the uniform bound on the exit rates in Lemma 3.5 


It turns out that the optimal size is controlled by the parameter (|2.12|i, and the crucial estimate in this 
context is ( 5.24| i for the probability that the trace process has changed well on the time scale 9l, where 
terms of the form and [3]^^ appear in bounds of the right-hand side. The best choice of (3^ in 
( |2.13| l leads to a bound in Lemma [33] which is small enough for the required estimates in the proofs of 
Propositions |2.1[[TT| and [33| as long as & > 20. The mechanism leading to this constraint is of a funda¬ 
mental nature, and it seems very hard to signihcantly improve this with the techniques used in this paper. 


The optimal choice of the second parameter ul in (2.121, ( 2.13| l controlling the fluctuations of the 
condensate is of a more technical nature, and arises from two conditions on the upper bound error of the 
capacity estimates in (7.12 1 . To obtain the upper bound, a test function is constructed that concentrates 
on wells and tubes (7.4 1 , which are subsets of A where the transition paths between wells are expected 
to concentrate on. The size of tubes and wells both increase with the parameter a, and increasing a 
therefore reduces the distances between wells and increases the corresponding capacity. At the same 
time, increasing tube size improves the approximation and therefore decreases the upper bound on 
those capacities. The optimal choice of results from estimates in Lemmas 7.4 and 7.5 and only 
leads to the weak requirement 6 > 5. 
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Possible extensions. We focus on symmetric, nearest-neighbour probabilities p{x,y) = \5y^x-i + 
^Sy^x+i in one dimension with periodic boundary conditions, two properties of which are essential for 
our proof. 

• Symmetry; this leads to reversible dynamics, which is a necessary condition for the potential the¬ 
oretic estimates on transition rates we use in Section]^ There is significant recent research interest 
on extended Dirichlet principles for non-reversible systems which involve double variational for¬ 
mulas Ha [TOl ED, see also references in CSl, Chapter 7. This has been applied to the totally 
asymmetric zero-range process on a fixed lattice in llJTl , but since a result in the thermodynamic 
limit requires much better control on error bounds an extenstion to this case would be a signif¬ 
icant technical complication. While a (non-optimal) relaxation time estimate for non-reversible 
systems can probably be obtained, we also make critical use of symmetry of the jump rates in the 
regularization step in Section]^ which is not obviously adapted to asymmetric situations. 

• Translation invariance; the results in 00 apply to zero-range processes without this property, 
leading to spatially inhomogeneous limit dynamics on a fixed lattice which are directly related to 
the choice of p{x, y). We use translation invariance in our proof for equilibration and also in the 
regularization step in Sections l^and]^ Specific simple examples of non-translation invariance, 
such as alternating p{x, y) or isolated inhomogeneities, can be treated as a direct extension of our 
result on a case-by-case basis. However, it is not clear how to formulate a result in the generality 
covered in 00, where the first problem already may arise when defining the limiting process if 
the probabilities p{x, y) do not have good scaling properties. 

Whenever the p{x, y) are symmetric and translation invariant and admit a well defined scaling 
limit of the capacities of a single random walk on A analogously to ( |6.6| l, our results can be extended 
without much effort to lead to limit processes with stationary, independent increments. General finite 
range symmetric p{x, y) in one or higher dimensions which scale to Brownian motion should all 
give the same result related to the corresponding harmonic functions of a single walker, appropriately 
modified on the torus. Note that in three and more dimensions these functions have a constant scaling 
limit leading to uniform displacement of the condensate, with expected logarithmic corrections in two 
dimensions. Also if p{x, y) has range diverging with L with well defined scaling limits for capacities 
our result can be directly adapted, including for example uniform p{x^ y) which leads to uniform 
condensate dynamics on the limiting torus. Due to the special properties of one-dimensional diffusion 
the case covered here is already one of the most interesting. Note also that even for finite-range p{x, y) 
the condensate dynamics will be non-local in all dimensions, in contrast to an analogous result for 
inclusion processes Il34l . 


In addition to the jump rates g{n) (2.3 i considered here, there are various other choices that lead 
to condensing zero-range processes (see e.g. ED and references therein). A well studied example is 
for rates with asymptotic behaviour g{n) ^ 1 + b/n^ with A S (0,1) leading to stretched exponential 
tails for the stationary measure EDEl. The lighter tail of the measure increases the depth of the 
metastable wells and leads to free energy barriers that grow sublinearly in the system size L, which 
is much faster than the logarithmic growth for the present model corresponding to A = 1. There¬ 
fore, even though an actual generalization would require considerable work, we expect that all our 
techniques can be applied and some estimates, in particular the ones in Section]^ should even get easier. 


The depth of metastable wells increases with particle density, and with the parameter b which deter¬ 
mines the tail of the stationary measure. In contrast to the restrictive conditions on b, our proof is robust 
in this system parameter and does not require any additional constraints on the particle density except 
p > Pc- As long as the excess mass N — pcL ^ y^LlogX it is known that the zero-range process still 
exhibits the condensation phenomenon on the level of stationary measures 0. It would therefore be 
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interesting to investigate, but is beyond the scope of this paper, in how far our results can be applied also 
for subextensive excess mass. 

Related work. The results presented here and other metastability results for zero-range processes 
only concern the stationary dynamics of a single condensate, corresponding to the slowest time scale in 
the system. In particular, these results only apply when the process is started on one of the metastable 
wells. Starting the process from a uniform initial condition with supercritical or diverging density 
leads to a dynamic formation of the single condensate on a different, faster time scale. This has been 
discussed heuristically in 051 [32l . where it was found that after a rapid formation of several large 
clusters, they exchange particles in a coarsening process on the time scale which leads 

to the formation of a single condensate. Note that the location of the clusters on this time scale is 
hxed. Recently the hrst rigorous result in this context has been obtained on a hxed lattice in the limit 
N —> 00 in jS). The coarsening dynamics takes place outside the set £ and is therefore approached 
with entirely different methods than the ones in this paper. For the condensation phenomenon in the 
inclusion process the dynamics of the condensate and the coarsening process take place on the same 
time scale, and both have been rigorously understood in 041 . 


4 Equilibration dynamics in the wells 

To establish an upper bound for the thermalization time scale of the trace process in a metastable well 
we follow the procedure outlined in ||9l. We introduce a process restricted to the well (called reflected 
process in Q), which is reversible w.r.t. the restricted measure /r, and estimate the relaxation time of 
this process. A general result from a can then be applied to yield a simple proof of Proposition 


4.1 Proof of Proposition |3.3| 

For each well £^, x G A, the restricted process is defined by the generator 

-/(^)) for all ?7 e . (4.1) 

CG£- 

As before, r{rj, () denote the jump rates of the full zero-range process, and jumps outside the well £^ 
are suppressed. Note that this is not equal to the trace process on £^, which has additional rates at the 
boundary of £“. It is easy to see that this process is irreducible on £^ and that it is reversible w.r.t. the 
restricted measure 




(4.2) 


The following estimate on the relaxation and mixing times of the restricted process will be used to prove 
Proposition |3. 3 1 


Lemma 4.1 The relaxation time tind the e-mixing time 
£^ are independent of x and bounded by 


^(e) of the restricted process TA [4.1) on 


frei < and fniix(e) < CL^{\ + L ^ log j) , 


(4.3) 


for some constant C > 0, depending only on the fixed parameters p and b. 
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The proof of this Lemma uses path counting techniques to establish a Poincare inequality for the re¬ 
stricted process, which is independent of the rest of this section and is therefore postponed to Subsection 

IQ 

We will use the following L? estimate for the ergodic average of a function that has mean zero on 
all wells, starting from the stationary measure fjf = p.[.|£’] restricted to the wells. The proof of can be 
found in ||£|, Section 3.1. 


Lemma 4.2 For every function /:£’—>■ K which has vanishing mean = 0 for all a; S A, and for 

allt>0 


E,. 


f{ri^(s))ds 


xGA 


(4.4) 


where = f^ei A the relaxation time on of the restricted process, which does not depend on x in 
our case. 


By translation invariance we can simply focus on initial conditions in a chosen well in the fol¬ 
lowing. To prove Propositionwe have to show p.ll| i, i.e. prove that 


E^ 


= Er 


l-OLt 


(r^(z) - r^{r]{s),£^)) (/(V'l(?7(s)) + z/L) - /(V’l(??(s))) 


zeA 
2 9^0 


0 (4.5) 


as L — 00 for all Lipschitz functions / : T — K and j] G £'^. For the total jump rate of the trace 
process into another well z ^ 0 , we have used the obvious notation 

r£(r],£^)= for all 77 e . 

Cee- 


With the definition of r^(z) in (3.101 we have 

r^{z) = for all x,zG A, 

and since the function / o is constant on all wells, the function 

hfiv) ■= Y {f{fL{v) + z/L) - f{fL 

2/0 


(4.6) 


has mean zero under p.^ for all x G A, independently of /. On every well £^, we can estimate its second 
moment as 


- 2^0 


2/0 


<Cf\\YrHGn E^^W^C-Q 


25.^0 


2^0 


1 logL 
L® log^ L Ol 


(4.7) 


Here, C/ is the Lipschitz constant of /. The last inequality follows from ( |3.13| l in Lemma 3.5 and 
equation ( |6.7| l, derived independently from capacity estimates in Section |7] The latter implies that the 
espected total exit rate 


0 L 4" A, with cl —?► 0 as L —>■ 00 . 

2/0 


(4.8) 
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With Lemmas Lemma 4.2 and 4.1 this implies that 
rteL 




hf{ri{s))ds 


< CCf t0L 


1 


,4 -- 1 


log^ L Ol 

with initial condition under the stationary distribution . Now 


= CCft 


L\ogL 


(4.9) 


sup E^x 

x£A 


([ df{r]{s))ds) <^E^x (f hf{r]{s))ds 

Wo ^ J LWo 


= L^fie 


tSL ^ 2 

hf{v{s))ds 


and therefore 


sup E^x 

xGA 


teL 


hf{ri{s))ds 


< CCft 


logL 


(4.10) 


Finally, we use Lemma 3.5 and our estimate on the mixing time in Lemma 4.1 to get for e of order 1/0l 


Ex 


nte. 


hf{r]{s))ds 


< Ex, 


^^mix (^) 


hf{ri{s))ds 


- Ex 


< C'/ll (fmix(e) + ef6»L) + sup ( E,x- [ 

II ~ xga V Jo 

<CCf({L^ + t) / , + Vf/logL) -^0 as L —oo. 

log L ' 

This finishes the proof of Proposition |3^ 


ptOL 

I hf{ri{s))ds 

^mix (^) 

2 \ 

hfivis))ds 


(4.11) 

□ 


4.2 Proof of Lemma 14.11 

We will derive an upper bound on the relaxation time of the restricted process in a well by proving a 
Poincare inequality, due to translation invariance it is enough to focus on the well f We will use that 
the stationary measure jjP outside the condensate location is essentially given by a product measure at 
the critical density in the limit L —^ oo. With /J.[f °] = ^[E]/L we can write 

/N = M°[(%,??A\o)] = 

where we use the notation p = (po, 77 a\o) for rj G to indicate the condensate size tjq and the bulk 
configuration pA\o outside the condensate. For fixed particle number N (i.e. under the measure p), 
Va\o obviously uniquely determines po- To simplify notation we identify measures with their mass 
function, and to avoid possible confusion in this section we will indicate the dimension of the product 
measure i/ in each term with a superscript. 


We interpret the product measure ^ in (4.12 1 as the stationary measure of an auxiliary system 
of L — 1 independent birth death chains with birth rate 1 and death rate g{n). Corresponding to our 


definition of the metastable wells in (2.12 1 , we restrict the state space of the chains to X = {0,..., /3 l}. 


Each chain has therefore the modified stationary measure 


-y^ = i^^l.\X] and 


[^] 


for all n G X . 


(4.13) 


Note that the state space of the auxiliary chains X^ ^ contains in particular all bulk configurations ?7 a\o 
for rj G E^, and recall that the product measure ly is stationary for the zero-range process. 
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(4.14) 


Lemma 4.3 There exists a C > 0 such that for all f : ^ 

Vav(/) < . 


and all L large enough 


Here 


L-l 


L-1 


= 


{CgX^- 1 : (iv-^c.,c) e£0} = {ceX^-i:^Cx<Pci + aL} (4.15) 


21 = 1 
-L-l 


is the subset of configurations in which are compatible bulk configurations in the well £^. For 

p G £^ we use the obvious extension f{r]) := /(?7a\o) with a slight abuse of notation, and the claim 
holds in particular for all functions defined on £^. 

Proof. Since p^{p) — = infc M°((/ ~ c)^), using the notation f = f — we have with 

( |442l i 


v.v(/) < M/ ) = ^ E , JV] t (-)*\o) . 


(4.16) 


Following a result by Doney 1241, = ./V] — Liy^[N — [pcL]] (l + o(l)) as L, N -G oo and 

N/L ^ p > Pc, and therefore for L large enough there exists C > 0 such that 

v.v(/)<cg^ ^ 


< 


ly^N - [pcL] - [aL]] ^ i-ir ^ 


(4.17) 


where we have also used p\£] —>■ 1 and monotonicity of n v^[n]. Since v has power-law tails (2.4i 
the first ratio converges to 1, and with the notation (|4.15 1 we get for large enough L 


Var^o(/)<C ^ p-^[C]PmBoiO , 


which finishes the proof with the definition of z/ (4.13 i. 


(4.18) 

□ 


The Dirichlet form of a single birth-death chain is given by 

= W [gp){f{n - 1) - /(n)) + lx{n+l){f{n + l) - f{n)) 

n —0 
Pl-1 

= +1) -/H)^ 

n=0 


(4.19) 


where we have used v[n] g{n) = v[n — 1] for n > 1 (2.4i and g(0) = 0. The Dirichlet form for L — 1 
independent chains is therefore given by 


L-l 


da_i(/) = ^ ^ p^-i[c](/(n-/(c))^ix(c. + i), 

X=1 


(4.20) 


where we write for the configuration where a particle is added to ( at site x, Cf = Cz + ^z,x- 
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Lemma 4.4 For all L > 1 we have a Poincare inequality for L — 1 independent chains, i.e. for all 
f : ^ K 




(4.21) 


For functions / = /Iso concentrating on bulk configurations of 8^ we can make a stronger statement. 


var,.-.(/iso) ^'^■'[c](/(r)-/(c)) iBo(r) 


x=l (^^=0 

:= ■ 


(4.22) 


This implies that the relaxation time for the independent birth-death chains is bounded above by /3£. 
This might seem like a crude bound, but it can in fact be shown that the relaxation time scales like that 
of a symmetric random walk na, even though the chains are driven to the origin and have stationary 
measure So our upper bound is sharp in the scaling /?£, and prefactors are not important for us here. 
Note that the stronger statement for bulk configurations ensures that in all terms the function / is only 
evaluated on B^, which is important to avoid contributions from the boundary of in the final estimate 
by the Dirichlet form of the restricted process in Lemma 4.5 below. 

Proof. We use a standard Efron-Stein estimate [TSll to bound the variance for L — 1 iid random 
variables. Writing A' = A \ 0 for the bulk sites 1 to L — 1 this is given by 


L-l 

Varpi.-i(/) < ^ i^^"MC](/(C) - /x(Ca'\x))^ , 

where we may choose any measurable function f^ : 
can simply choose 

MCa'\x) :=/((Cav-0)) . 


(4.23) 


L To show the first statement (4.21 1 we 


(4.24) 


where as before (Ca'VxjO) denotes the configuration ( where is replaced by 0. Using the Cauchy- 
Schwarz inequality 


Cx-i 


(/(C) - /(Cav^o)) < C. E (f(CA'\x,i + 1) - /(Cav.O) , 

this leads to 


(4.25) 


1=0 


L-l 


Pl Cx -1 


Var,.-i{/)<E E EC-^MCx] E(/(Cav>^+1)-/(Cav-0) • (4.26) 


^A'\x 


Cx=o 1=0 
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Reordering the sum and using that n i— ni^^[n] is monotone decreasing, we get 


L-l 


/3i-l 


I3l 


Varpi.-i(/) < ^ '^[CA'\x]'^{f{CA'\x,i+^)- f{CA>\x,l)) CxI^^iCx] 


‘•A'\x 


L-l 


1=0 


Pl-1 


Cx-1+1 


E E(/(Cav>^+1)-/(Cav^0) 


a;=l «A'\x 


i=0 


</3i/4 






(4.27) 


The same argument works when we restrict to functions /Iso. Note that G implies 
(Ca'\:c ) 0) G i3°, and also all configurations ap pearin g in the Cauchy-Sc hwarz decomposition in 
(4.25 I are in _B°. Therefore, restricting the sum in (4.23 i to G leads to (4.22 1 with a completely 
analogous computation, which finishes the proof. □ 


To finish the proof of Lemma |4~T| we will bound the Dirichlet form of L — 1 independent birth-death 
chains restricted to by the Dirichlet form of the restricted process using a standard path counting 
argument. The bounds we get here are certainly not optimal, and one of the reasons for our conditions 
on the parameter b. 


Lemma 4.5 There exists C > 0 such that for all f : 


. (cf. {4.15\) and all L large enough 


<.(/) < 


where B^^f) is the Dirichlet form of the restricted process 


(4.28) 


(4.29) 




Proof. For each / : — >■ M we write /(p) := /(77 a\o) for its unique extension to 77 G as before. 

The Dirichlet form of the restricted process is simply given by 

^°(/) = T E E ^ 

TjeS° xeA 

+ (/(^-.-i) _ f(r,))heo{r,^’^-^)) , (4.30) 

since all jumps leading outside are suppressed. We change the summation to the set 

£%_i = {r] € ■. rj^ £ for some z G A} , (4.31) 

and use nV l^]9{'dx + 1) = Ml n-iIv] "'i’-f* '■f*® canonical measure for iV — 1 particles, which follows 
from ( |2.4| l. The Dirichlet form can then be written as 

'?6£Sr-l ^eA 
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With the above notation the restricted Dirichlet form ( |4.22[ ) of the independent chains can be written as 

r]eS° xeA 

= E E - ( 4 - 33 ) 


where we used the same change of summation variable as above in the second line, and the fact that 
f'^“^[(?7°)A\o] = We decompose the transport of a particle from the condensate site 0 to a; 

into nearest neighbour jumps and use the Cauchy-Schwarz inequality on the telescoping sum to get 




(4.34) 


y=o 


for all X, since the longest path of a particle is clearly bounded by L. We can bound every such term 
in (4.33 I that way, and as long as rjy + 1 < for all y = 1,... a: — 1, all ry*' G and the terms in 
the sum (4.34 1 correspond to ’allowed’ transitions that appear also in I?°(/). The ’flow’ of an allowed 


transition is the number of times it appears in (4.33 i using (4.34 1 , and summing over all target positions 
X this is bounded by L. 


On the other hand, if there exist sites 0 < yi < ... < < x, m > 0, with r/y. = /3 a, the generic 


path in (4.33 1 contains non-allowed transitions and has to be re-routed, increasing the flow of certain 


allowed transitions. To bound this increase, we introduce the notation 


cr^ry := ty + for y, z G A , 


(4.35) 


where Sy is the configuration with a single particle at y. The path corresponding to the sum (4.34i can 
then be represented by the equation 

ry“ = cr;^_i •• •cr?cro^° . 

If there is an isolated site y with rjy = /3a and ryy±i < 13 re-route the path of a particle from y — 1 
to ?/ -f 1 from 


ry^+i = al+^al_^Tf ^ 


to = ay_iay'^^r]y ^ 


Instead of moving the particle to site y, it remains in site yy — 1 and a particle moves from yloy + 1 first. 
In the next step the particle follows from site y — Ito y reaching 7y^+^ only via allowed transitions. If 
there is a block of k consecutive sites with rjy = ... = rjy+k-i = Pl, the re-routed path of a particle 
from y — Ito y + k along valid transitions is 


„V+k _ V ... y+k-l y+k y-l 
V - ^y-l <^y+k-2'^y+k-l^ 


(4.36) 


Possibly combining re-routing over several blocks of sites with occupation number Pl, we associate a 
unique particle path from 0 to a; to each base configuration p G Xa.at-i, using only allowed transitions. 
The flow of a transition is then multiplied by the number of associated base configurations that use it 
for some x G A, and every transition with multiplicity higher than 1 involves at least one site with 
occupation number /3a. Denote by (j — (j' one of the transitions along the path in (4.361, then one 
associated base configuration is obviously ry, and another one is given by the minimal configuration 


C A C' := (C. A C : ^ G A) . 
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Note that for all transitions C ^ C along the path in (4.36 1 we have for the maximal configuration 
c V C' := (C. V c; : ^ G A) = = r, + <5,_i + Sy+k ■ 


It is easy to see that any possible base configuration associated to a transition ( —> CAn (4.36 1 has to be 
of the form 

jjV-^.v+k _ for some z,z' £ {y — 1,... ,k + 1} . 

In many cases not all of those base configurations contribute (or are even in £), but this p rovides an 
upper bound of (/c + 1)^ for the flow multiplicity of transitions C ^ C along the path (4.36|l. 


For any base configuration p summed over in (4.33 i there are at most of order L//3i, sites with oc¬ 
cupation number /3l, and therefore, the multiplicity of the flow along any a llowed transition is bounded 
by C{L//3l)‘^. Together with the generic flow bound of order L and (4.34 1 this implies that 


Using again ( |4.12[ ) and the same approach as in the proof of Lemma [43] we can bound 


,L-1 


-L-1 


[^A\' 


u^[Sl=N- 1] yL.N-i[£] 


[^^AXo] L iy^[r]o] 

-c ^ uO r„i 


4iv] 


(4.37) 


for a suitable C > 0. We have used again that v^[Sl = A^] = Lv^[N — [pcL\\ (l + o(l)) 
monotonicity of n i—)■ i/^ln] and pL^pf-i[£] —>■ 1. Since 


A[N] 


we end up with 






which finishes the proof of Lemma [43] 


□ 


Together with Lemmas |4.3| and |4i4| we get a Poincare inequality for the restricted process, i.e. for all 
f : £° ^M. there exists C > 0 such that 


Var^o(/) < CL^V^if) . 


(4.38) 


Therefore, the relaxation time for the restricted process is bounded by frei < on each well £^, and 
by a standard result lf3^ this implies for the e-mixing time that 

(4.39) 

□ 


fmix(e) < -Lei log (e miii/x°[77]) < CL^{l + L ^ log(l/e)) , 
since pP[r]] is at most exponentially small in L. This finishes the proof of Lemma |4.l| 
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5 Uniform bounds on exit rates via coupling 


To derive the uniform bound of Lemma 3.5 on the exit rate out of a well for the trace process, we will 
construct a coupling of the zero-range process with a growing number of birth-death chains. The number 
of chains increases only linearly in time, and the coupling ensures that in the event of changing well, 
at least one of the chains has grown a condensate, the probability of which can be controlled directly 
from metastability results on a fixed size lattice in jT]. We will need additional control on how much 
time the full process spends outside the well which we achieve by using mixing estimates on larger 
wells containing the original ones. This is derived first in the last subsection, together with a proof of 
Proposition |2. 1 1 on substitution by the trace process. 


5.1 Construction of the coupling 


We construct the coupling that will be used in the next subsection to prove uniform bounds on exit rates 
from wells. Let (C(f) : < > 0) be a birth-death chain on the state space Jf = {0,1,...} with birth/arrival 
rate 1 and death/departure rate g{Q as given in characterized by the generator 


CfiO = (/(C +1) - /(C)) + 5(0 (/(C - 1) - /(C)) • 


(5.1) 


Note that the boundary condition at = 0 is included with p(0) = 0, and this chain has stationary 
measure i/ as given in (2.4 1 . For some fixed e G {0, p — pc) denote by 


VL ■■= {p- Pc- e)L and Ty^ := inf{f > 0 ; ({t) > yL} 


(5.2) 


a size-dependent level that has to be crossed to grow a condensate, and the associated hitting time. 
Lower bounds for the hitting time Ty^ are typically of order and can be derived by direct 

computation. 


Lemma 5.1 There exist constants Ci,C 2 > 0, such that for all initial conditions (^(0) = Co G 
{0,1,..., B]f\ with 1 <C Bl L, we have that 


EcATyc] > CiOl 


b-1 


and [Ty^ < t] < C 2 


(5.3) 


Proof. It is easy to show that the expected hitting time Ty = [Ty] with a; < y € N of a birth-death 
chain with birth rates h{(), death rates p(C). and stationary measure u is given by 


y-i 

r;=E,[T,]=^ 


C 


„=o 


(5.4) 


For a reference see e.g. We have h{Q = 1 and due to monotonicity of u in this case we can use 
simple integral bounds for sums. We get 


1 fC 

u[n] > — / u~^ du = 

n=o -^1 


1-C 


l-b 


Zc{b- 1) ’ 


which then analogously leads to 


>E —> 


C=x 


6-1 


{y - 1)''+^ - 
(h-l){b + l) 


y"^ — 

2(&-l) ■ 


(5.5) 


(5.6) 


22 










With x<BL<^y = yL^L this directly implies the first statement. 


To derive the second statement, for a given Co we couple the chain with a modified chain C^( ) that 
cannot jump below Coj i-C- it has death rates 


g'iO = 5(C) for c > Co and ^'(Co) = 0 . 


It is clear that the chain ('{■) will reach y^ before the original one, so its hitting time Ty^ will provide 
a lower bound for . Furthermore, the lowest hitting time is clearly achieved for Co = and we 
can focus on this case. Since the point Co = is the left end of the state space for the C^(-) chain, the 
Markov property implies the following sub-multiplicity, 

> s] < . 

Integrating over s and re-arranging yields 


P),< f] < f/E),JT;J , 

and it remains to estimate the expectation of T'y^ from below. Note that the chain C^(') has the same 
stationary measure v restricted to C > Bl, changing the normalization to z'. The latter cancels in (5.4 1 , 
and we simply have to adapt as 


Y . ^ 77 


r-C _ ^1-fc 

u-^du= ^ ^ 


'l=Br 


JBl 




Analogously to (5.6 1 this implies 


VL-l ^ 


C=Bi. 


6-1 


'7> 


for a suitable constant C 2 > 0, finishing the proof. 


□ 


Now let us fix a configuration rj G For the original process (p(f) ■ t > 0) the occupation 
numbers outside the condensate rjxit), x = {1,..., L — 1} are birth-death processes with Markovian 
departure processes at rate g{r]x), but non-Markovian arrival processes that depend on the neighbouring 
occupation numbers for each site. Conditioned on the configuration ry(f) at time t, the arrival rate at site 
X is given by 

ax{r]{t)) := ^g{r]x-i{t)) + ^ g{r]x+iit)), (5.7) 

and if both neighbouring sites are occupied this rate can be as large as 2^. In order to dominate gx{t) by 
a Markovian birth-death chain with arrival rate 1 to apply Lemma [5T| we couple it with an increasing 
number of chains. At any given time, at least one of those chains will dominate rjxit) and lead to an 
estimate for the probability of leaving the well 


Let TO > 2^ be the smallest integer greater than or equal to the maximal jump rate g{2) = 2^ (2.3 1 
of the zero range process. The coupling as described below is applied for all times t > 0 and sites 
X = 1,..., L — 1, and is illustrated in Figure[T]for the simplest case to = 2. To each site x we associate 
an infinite number of birth-death chains (C^ (f) ■ t > 0) with generator (5.1 1 , where we index the chains 
by vectors of variable length of the form k = (fci,..., fc„) with ki € to} and n = 1 , 2 ,.... 
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This corresponds to indexing the chains by the nodes of an m-ary regular tree TZm without root, with 
generations indexed by n, and we write k G 'R-m- At any given time, each chain C,'^(t) in the tree 7j is 
an identical copy of its unique parent chain, 2) evolves independently of T]x{t) and all other chains, or 
51 is associated to rix{t) as described below. The assignment of each chain to one of these three groups 
changes in time, and we denote by 

Cx{t) := {k G TZm '■ Cx (0 associated to r]x{t), and no ancestor of Cx (0 associated to ?7x(i)} 
1x{t) := {k G TZm ■ Cx (0 evolves independently} (5.8) 

the index sets of chains which are not identical copies of their parent. At any time t > 0, the number of 
chains in Cx{t) is |Ca;(f)| = m, for all sites x. 


Initially, we set the m chains in generation n = 1 equal to rjx, i.e. Cx{Q) = {(1), ■ ■ •, (w)}, = 0 

and all other chains are identical copies of their parent. We use identical initial conditions, that is, for 
each site x 

= ... = C}™)(0) = px(0) G {0,1,... ,/3i} , (5.9) 

and our coupling will ensure that rix{t) < C^{t) for all k G Cx{t). For the departure process of 
associated chains with k G Cx{t) we simply use a basic coupling for all x and t > 0, i.e. particles 
in Cx(0 leave together with particles in ryx(f) with probability g{Cxi^))/ 9 {Vx{t)) < 1 for r]x{t) > 1, 
and they additionally leave, independently of particles in r]x{t), at rate g{Cx{t)) ~ g{Vx{t)) in case this 
quantity is positive for r]x{t) < 1- Note that the departure dynamics preserve the order 

ilx{t) A Cx (^) f°r all k G Cx{t) and t> 0 , (5.10) 


and we will couple the arrival processes in such a way that this is true also for the full process. To 
achieve this, we change the structure summarized by the sets Cx{t) and Ix (t) at every jump event on the 
arrival site. When a particle arrives at site x in the p-process at time t we pick one of the m chains in 
Cx (f) uniformly at random, add a particle to all of its m children (which up to this point have evolved 
as identical copies of their parent), and disassociate the chains in Cx{t—) so that from this time on they 
run independently of gx{s), s > t. That is, sample k* ~ {7 {Cx{t)), and let 


Cx(f) = {(k*,l),(k*,2),...,(k*,m)} , Ix{t)=Ix{t-)UCx{t-) . 

So far the coupling leads to an effective arrival rate of ax{r]{t))/m < 1, ax{ri{t)) as in ( |5.7| i, for each 
associated chain, which is typically strictly smaller than 1. To each associated process we independently 
add particles at rate 1 — Ux ivit))/m, leading to a total arrival rate of 1 as required. Therefore all chains 
(Cx (f) : f > 0), k G TZm, have the desired marginal dynamics of a birth-death chain with generator 
(IH)- Note that the total number of particles in the associated chains is not conserved and is growing in 
time, but the main point is that the coupling fulfills (5.101. 


This coupling construction leads to increasing sets Ix (t) of independently evolving chains, but at 
any time there is only a finite number of chains which are not identical copies of their parent. This 
number grows only linearly in time with high probability, and we will use this in the next subsection to 
prove a uniform bound on the exit rate from a well. 


5.2 Proof of Lemma 13.51 

For the trace process (? 7 ^(s) : s > 0), ^7^(0) = 77 ( 0 ) G to reach another well before time t > 0, one 
of the occupation numbers for x 0 has to grow larger than = {p — Pc — s)!, e G {0,p — pc) given 
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Figure 1: Coupling procedure illustrated for m = 2, associated chains are shown in red and are encircled, 
independent chains are shown in blue in a rectangular box. Arrows indicate identical copies. Initially (top) only 
the chains in generation n = 1 are associated. After the first particle arrives at site x (middle) both chains in the 
first generation become independent, and the m descendants of the the chain get associated. This process is 
repeated after the second particle arrives (bottom). 
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in (5.21, before the full process (p(s) : s > 0) has spent a total of time t in the well f By construction 
of the above coupling, this implies that at least one of the associated chains with k S TZm and 
must also have reached this level. Again, it is important to note that for each x 0 and s > 0 there is 
only the finite number of chains in Cx(s) U Ix{s) C TZm which are not identical copies of their parent. 
Since all the chains associated to site x are birth-death chains with generator ( |5.1| l, we can use Lemma 
|5.1| to estimate the time it takes for the chains to reach level y^. 

Note that the dynamics of the birth-death chain {C,{t),t > 0) does not depend on the parameter L, 


hence Lemma 5.1 in fact says that for any parameter K ^ oo, letting Tk ■= inf{t > 0 : C(f) > K} 


and for an initial condition < Bx, 1 ^ Bx K we have 


Co 


[Tx <t]<C 


tB 


6-1 

K 


Rl+b ■ 


(5.11) 


We apply this bound to gain control on the time it takes the process {rj{s) : s > 0) to exit a larger 
well 8^ D 8^, which we choose as 


8° := [rj e Xl,n ■ Vo > ^ - PcL - ul, Vy < (j^L for all y ^ O}, with (pL ■= 


L 


logL ■ 


(5.12) 


Clearly, the larger wells dehne an analogous partition of the state space that uniquely characterises the 
condensate location for sufficiently large L, so t he de finition of the restricted process ( |4.1| l in Section 
l^can be adapted to the extended wells. Lemma |4 a] directly applies to 8^ replacing 4>l in the 

derivation of all estimates, and therefore yields the following bounds on the relaxation time tre\ and the 
(5-mixing time finix(^) of the restricted process in 


frei < CL‘^ and ts := fmix(^) < CL^ {L + log(l/5)). 


(5.13) 


Let us denote by T the exit time of the process (p(s) : s > 0) from the set 8^ and for p G let us 
write 


B{r]) = maxrjx 

x^O 


Lemma 5.2 There exist constants C, 7 > 0, such that for any 7 € £’0 ond t > 0 

Mf<t] ' +Le-^*. (5.14) 

•Pl 

Proof. On account of the coupling argument of Subsection |5.1| if we consider the system of birth-death 
chains (C;^(s) : s > 0,a; G A \ {0},k G TZm) with initial condition associated to 77 G 8^, so that in 
particular Cx (0) < B{r]), x ^ 0, ]<. G Cx{0), then 

{T < f} C |3k G TZm, X ^ 0 and s G [0,f] such that C7s) > 


so 

^r,[T < f] < PcO) for some s G [0,f], x ^ 0, k G C 2 ;(s) U Ja:(s)]. 

The number of chains in Cx {t)Ulx (t) increases with arrival processes of particles which are independent 
Poisson with bounded rates, so the probability that |Ca;(f)| + |Ta:(f)| > Ct decays exponentially in t for 
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C large enough. Applying ( |5.11| l to each of these chains independently with K = (pL we obtain 
IP 77 [T <t] < + IP 77 [|Ca; (f) I + \Tx (f ) I > Ct for some a; ^ O] 

<C +Le , 


(5.15) 

□ 


which completes the proof of the Lemma. 

We turn now to the proof (|3.13|l in the statement of Lemma 3.5 Let us first note that for any f > 0 


Pr; 0 for some s G [ 0 , f]] > ^1 — e j 


E 








(5.16) 


where the ratio on the right is the probability that the first jump out of the configuration p is to a 
configuration in the complement of Let now f = j. Since 

XI ^ diVx) < mL 

^^ 7 ] xGA 


we get that 


X < CLF^[tjjLiv^is)) ^ 0 for some s < -], 


(5.17) 


where 


C = 


- - = sup --. 

^ ^ a:G(0,m] ^ 


We next prove that there exists a constant C > 0 such that 
sup [tpLiil^(s)) ^ 0 for some s < y] < C 

Lj 


ri&E° 


L6 log^L’ 


(5.18) 


which together with ( |5.17| i immediately implies the assertion of Lemma |3.5| For this, write 
{ip l{t 1^ {s)) ^ 0 for some s < ^} = {r]{s) G £\£^ for some s G [ 0 , S'l-i]} 

= { j lfo(77(s))ds < -}, 


(5.19) 
and 


2.2 


where St := sup{s > 0 : 7^ < f} is the inverse of the local time Ts in 8^ defined in Section 
T£\£o is the hitting time of the process (77(3) : s > 0 ) on £ \ To control the probability 

of the event in ( 5.19[ ) we define an intermediate, deterministic time II which is small enough for the 
process rj{-) started at 77 G £° to remain in £° by with high probability, but large enough for the 
restricted process on £°, which we denote by (^(s) : s > 0), to have mixed. Precisely, we require 


— > CL'^ for some k > 0 and 
ts 


sup P,, [f < fi] < 


C 


as L 


n(s£« ■ - log L 

By (5.13 I and Lemma [0| these conditions are simultaneously satisfied if is chosen as 

y(b-10)/2 

(logL)^ 


(5.20) 


(5.21) 
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In fact, the condition b > 20 arises from the need that ti, ^ ts in (|5.20|l. Define now 


A = { [ l£o {r]{s))ds < y} 2 {'0 l(^^(s)) ^ 0 for some s< 
Jo ^ 


For r] G we have 

Pr, [-4] < p^f < fi] + n {T > fi}] 


<Frf[f < ti] +P^[^ l£o(7?(s))(is < ^; T > fi] 

< P^ [f < fi] + Pr,(ii) [A]-, f >tL , 


(5.22) 


(5.23) 


where the third line follows from the Markov property. We may couple a restricted process to £^, which 
we will denote by (^(f) : f > 0), to the original process, so that they jump together up to time T. If 
P^™^) is the coupling measure with marginals P^ (original process) and P^ (restricted process), then 


E. 


- rj{t 


) [yf]; f > tL 


= E^^ 

[-4]; T > tL 

= E““P 

{v,v) 

< Er") 


= 

P 



T >iL 


We now substitute this estimate in ( 5.23| to get that for any rj G 


Pr; [-4] < P,; [T < fi] + E^ 


■^Ol) 




< P^ [T < fi] + P/iO [^] + dL, 


(5.24) 


where d/, is the total variation distance between the distribution of ^(t^) and the invariant measure pP 
of the process (^(s) : s > 0). By & pP is simply the invariant distribution /r of rjp) restricted to 
that is 


A°[e] = l£o(0 


p[£0] 


and by (4.35) in lf39l , choosing d < 1/2 we have 


dL < (2(5)'^'. 


When r] G ( |5.20[) ap plies and, in view of (5.22 1 , (5.24i and the observation following (5.18 1 , the 
assertion of Lemma [33] reduces to showing that 


Pa« [^] < 


c 

L6 log^L' 


(5.25) 


To estimate P^o [,A] we partition as y £ 0,2 y £0.3^ where 


£04 = {^ g £0 . > iv - + ^, S(p) <^}C £\ 

£04 = {^ g £0 . Bi-n) < xl} \ and = {,y g fO . B{n) > xl}, 
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for some xl to be determined later. We may now write 


sup P4^] + sup (5.26) 

r)G£0-i j7G£°-2 

and use ( |5.24| l to estimate the middle term of the sum on the right hand side. This gives 

{l-fL°[£°’^])Ff,o[A]< sup P4^]+( sup F^[f<tL]+dL)fi°[£°’^]+fi°[£°’^]. (5.27) 

r,e£OA Tie£°’^ 

By Proposition |9^ we have 

il°[£°’^]<CLPl-^ and il°[£°’^] < CLxl-^ 
and by Lemma [5^ 

sup P4T<f^] 

r,PF0,2 (hr 


When rj G £^’^ for the process to leave at least particles need to be moved, either away from 
the condensate or onto the site that already contains ^5^/2 particles. If fewer jumps occur by time P, the 
process will necessarily remain in £^ up to time and the event A will not be realised. Since the total 
rate at which jumps occur is bounded by mL, the probability that at least /3i/2 jumps occur by time ^ 
is dominated by P[2f > where X has Poisson distribution with mean m. Hence, 


sup P,, [yf] 
^^£■0,1 


< 


E 

fe>/3i,/2 


fc! 


< 


m 


/3i/2 


{/3l/2)!’ 


which decays faster than any power of L If we put together the preceding estimates, (5.27 1 gives 


Ff,o [yf] < CL 


Lti/3 


2 ol-h 




The estimate (5.251 now follows, if we optimise the preceding bound by choosing 


XL = 


1 + fcflfe-l \ 2(6-1) 

L 1 




Ltl 


= {LAog^LY'^ 


and this concludes the proof of Lemma [33] 


□ 


5.3 Proof of Proposition |2jJ- Replacement by the trace process 

We will make use of the extended well f ° and the exit time T of the process p(-) from as well as 
the process restricted to denoted by ^(-j, and its 5-mixing time ts- All are defined in the preceding 
subsection. We begin the proof of Proposition |2. 1 [with the following estimate. 


Lemma 5.3 For any u > 0 we have 

(ts+U 


sup 

r)e£° 


1a(?7(s)) ds; T > 4 < u ( y + ^ -i^[A] ) 


1 - 


(5.28) 
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Proof. For ry S ^ let us define W{ri) = 
property we have 

rts+u 


'-./O 


1a(?7(s)) ds 


. For any rj G using the Markov 


Er 




lA{r]{s)) ds; T > ts = Eri[W[r]{ts)); T > ts] < E,^[W[^{ts))] = Wdn^, 


where tt^ is the distribution of the random variable ^{ts) under P„. The inequality above follows by 
coupling ^(•) to ? 7 ( ) up to time T, just as in the argument following (5.23 i. Then, 

f WdiTr) < S sup fF(C) + f W dpP = 5 sup W{^) H- f W dfi 

J fpro J fefo u.\£°\ Jeo 


= 6 sup kF(^) 
ie£° 


f[U 


Jgfo fj.[£° 

W dpL 


xeA ' 




<i5mH-/ Wdfi = u(6-\ - ^^ J , ) 


i-m[a] Jx^ 


1 - n[A] / 


The first line follows from the definition of ts and ( |4.2| ), the second line follows from translation 
invariance of the dynamics and the set A and the final step uses the invariance of p under the dynamics 
ofr7(-). □ 

It is now straightforward to prove the replacement by the trace process. 

Proof of Proposition|2.H In view of Lemma [53] we have 


Er 


/•ter 


lA(T?(s))ds 


= E^ 


rtdr 


'-Jo 


1a(??(s)) ds; f <tL 


■ E.,. 


ntOT 


^Jo 


1a(?7(s)) ds; f > II 


<tOLEr^[f<tL\+ts+Er^ / 1A (??(s)) cfs; T > 


i-te. 


‘-Jts 


<ts+ t0L 




+ S + 


f"[A] 

1 - /r[A] J 


Dividing hy 9 l, the assertion now follows from ( 5.20| ) and Corollary |9.5| which provides a vanishing 
bound on p,[A]. □ 


6 Regularization and inter wells dynamics 


In order to get matching upper and lower bounds on the transition rates of the auxiliary process ( |3.9| l, 
we have to replace it with a regularized version on a renormalized latt ice, making use of Lipschitz 
continuity of the test functions /. In fact, we will show Proposition 3.4 for all / G C^(T,R), which 
can be used to uniformly approximate Lipschitz functions. 


Sections [^ and [^ are independent of the rest of the article, the proofs here rely on some invariant 
measure estimates provided in Section]^ The only restrictions arising from the results of these sections 
on the parameter b are given by equations \6A\ and ( |6.13| l, resulting on a lower bound 5 > 5. Clearly 
this is much better than the condition & > 20 following from the uniform bounds on the exit rates 
(Lemma [3.5| l, the estimates required to prove equilibration in the wells (Proposition |3.3| l and tightness 
of the condensate dynamics (Proposition [XT]). 
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For future use where a different scaling may be needed, we choose to keep ul and and other 
quantities derived from them explicit, instead of replacing them by the values obtained with the choices 
in( |2l3| ). 

6.1 Rate estimates from capacity bounds 


Generalizing the rates r^[z) given in (3.10i to non-empty subsets Ai,A 2 C A with Ai Cl A 2 = 0, we 
write 


(Ai, A 2 ) = -rj-r r^{y - x) , 


1^ 


( 6 . 1 ) 


y^A2 


using that under fi, conditioned on p G = UxgAi location of the condensate is uniformly 

distributed in Ai. In this notation we can identify r^(z) = r^({0}, {z}). Following ||6l. Lemma 6.8, 
we have 

r^(Ai,A2) = 

^(^ca.p^{£^\£\£^^)+cap^(£^^,£\£^^) - ,£ . (6.2) 

Here cap^ [£^,£^) denotes the capacity of the trace process between the sets of wells £^ and £^ for 
any A, B C A with A n B = 0. By Lemma 6.9 in ||6l, the trace process capacities satisfy 

y[£] cap^(£■^£:^) = cap(£’^,£:^), 

where the latter are the full zero-range process capacities. By Corollary 9.5 in Section 
m[A] < CL(^a\ ^ + j3\ ^), hence we may replace and cap^(-,-) throughout by the full 

zero-range process invariant measure and capacities /i[-] and cap(-,-), at the cost of an o(l/L) 
multiplicative error. 

Propositions |7 .1 1 and |7 .2| provide the following lower and upper bounds on capacities between com¬ 
plementary sets. For each A C A we have 

6»Lcap(£:"^,£’\£:^)>iJ(l-e^) ^capA(a;,y) , 

y^A 

9Lcap{£^,£\£^)<H{l + eL) capf,{x,y) , (6.3) 

xeA 

y^A 

with a constant H = H{b, p) = )i+bz ^^d error terms 


eL = C{LPl-^+ L-^+Le and el » > 


|A|L''+3 


A’ -af-^{L-\A\y 


(6.4) 


With our choice of aL and Pl in (2.131 ei, ^ ej^, and using (6.2 1 this leads to 


0Ly[£^ 


^ r‘^{Ai,A2)< capA(x,y)-f eif ^ capf^{x,y) + ^ cap/^{x,y) 

x^Ai ' x^Ai x^A^ 

y£A2 y^Ai y^A2 

Ai 1 / 

^ r'^(Ai,A2)> capA(x,y) - eL( capA(x,y)+ capA(x,2/) 


xGAi 

y^A2 


xGAi 


x€A2 

y^A2 


(6.5) 
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The bounds are given in terms of capacities of a simple symmetric random walk on the lattice A, which 
are 


capA(x, j/) =capA(0,j/ - x) = l/df^{Q,y- x) 

1 1 
\y-x\ L-\y-x\ 

using the standard embedding {0, ... ,L — 1} of the torus A in Z. In the following we will use a 
double embedding {0,..., L — 1} and {—L, ..., —1} shifted by —L, where we identify sites 0 and 
—L. The combined use leads to a more intuitive notation of the regularization procedure which is 
formulated symmetrically around 0, and should not cause any confusion. Note also that we write \x — 
y\ G {0,..., L — 1} when x, y are chosen consistently in the same embedding, and that this is not the 
distance c^a on the discrete torus A. Some particular bounds we will make use of later are 



\y - x\{L - \y - x\) L 


9Lr^{z) < Hl(^ 


1 


L-1 


z\{L- \z\) 


+ 2eL 




xLJ 


(r 


and 6l r^(z) < C(l + e^) log A 
2^0 


1 



(6.7) 


using that = (1 — p,[A])/L = 1/A(1 + o(l)). Note that 
from the first; instead, we exploit the fact that = 

directly. 


the second line in ( |6.7| ) 
r^(0, A \ 0) and apply 


does not follow 
Proposition 


7.2 


6.2 Proof of Proposition |3.4| 

To establish p.l2|i in Propositionwe will show that 


sup 

xGA 


OLC^fix/L) - C^f{x/L) 


0 as L —>■ 00 . 


( 6 . 8 ) 


To achieve this we have to regularize the auxiliary generator (3.91, since using bounds of type (6.7 1 
directly will lead to diverging error terms. Fix an intermediate scale i depending on L, which is a 
sequence of integer numbers such that 


^—>■00 and i/L —>• 0 as L —>• oo . (6.9) 

We partition the lattice A into subsets (or boxes) of size 2i + 1, with only the box Vq centred in the 
origin 0 being of larger size. For this we fix a second sequence £ such that 

f/f —>■ oo , i/L^O as L —oo , (6.10) 

and take Vq = {~i^ ...,£}. We choose £ such that the remaining lattice A \ (/q can be partitioned into 
boxes of size 2£ + 1, i.e. there exists a sequence of integers M such that 

L-{2£ + l) = M{2£ + 1) . (6.11) 

The integer sequence M is bounded above by the real-valued sequence M := L/{2£ + 1), which 
characterizes the asymptotic number of boxes, since M/M —1 as L —>^ oo. The partition of the lattice 
is then given by 1^ and the boxes 

V^:={yG A:\y-£-(2£+l)m\<£} , m=l,...,M. (6.12) 
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A choice of scales consistent with ( |6.9| l and ( |6.10| l which is sufficient for the proof is 

(. = 0{aL^o^ L) and £ = 0(ai log"* L) as L —>• oo . (6.13) 

This is possible with our choice of a/, in (|2.13|l, compatible with ( |6.4| i, and implies in particular that 


eLATlogi—T'O and — logL—^0 as L ^ oo . 

Lj 


(6.14) 


Here we have used that we can choose the error in (6.4i as ^ log A, since the size of the 

sets of wells we consider is |A| < f. In fact, a quick computation shows that among the sequences 
Ul = 7 > 0, and £ satisfying ( |6.9| l, ( |6.10| l, the choices ul = Ls + se and I = ul log^ L provide 

the smallest possible (up to powers of logL) which is consistent with ( |7.12[ ), and require a lower 
bound on the parameter b of only b > 5. 


To regularize the rates (3.10 1 , we rewrite the generator of the auxiliary process (3.9 1 as 
= -/(!)) 


zeA 

1 




2 £+ 1 


yeVo zGA 


L 


2i 


y&Vo zGVb 


M 


'EE E 

m=i yeVo zeVm 


r^iz- 


!/)(/( 


x + z - y 




(6.15) 


=-.Ctfi^/L)+C^f{x/L). 

In the hrst step we used translation invariance of the rates and averaged them over the box Vq = 
{—£, ...,£} of size 2£ + 1 around the origin. For the hrst term we use the Lipschitz property of / 
with constant C/ to get 

6»L/:^/(a;/L)<C'/—^ ^ ^ ^ r^{z) 

zeVoV&Vo O/zGVb 


< <^^(1 + eL)logT , 


(6.16) 


where we used (6.7 1 in the last line. So with (6.141, |£^/{a;/L)| —0 as L —>■ oo uniformly for x G A. 
The second term of (6.15 1 can be split as C^fix/L) = (xjV) + C^sfix/L) with 

'x + i+ i2£ + l)m' 


M 

C-^Afix/V) := ^ ^ ^ r^(z -?/)(/( 


X + z — y 


2 £+ 1 

m=l yeVo ZeVrr, 

M 




7e/Wi) := E ’■‘(ro. r„)(/( -^ + '' + f+ *>’" ) - /(f)). 

m —1 

For the first term we can use the Lipschitz continuity of / with constant C/ to get 


M 


9/ _|_ 1 

\cLf{x/L)\ < Cf Y. . 


(6.17) 


(6.18) 


m —1 
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Using (6.5 I for complementary sets Ai = Vq and A 2 = A \ Vb we get 




m—l 


xGVo 

J/^Vq 


< CK(1 + sl) ^ < C log L , 


(6.19) 


which leads to 


&Lj£Lf(^/L)l<Cf(2£+l) 


log A 


0 as A —>■ (X) uniformly in a; G A . 


( 6 . 20 ) 


In the main contribution £^d in (6.17 1 we can use (6.5 1 to obtain for the renormalized rates 


dLr^{Vo,Vm) = ^( ^ ^ capA(0,2;-?/) + eL7^L(Uo,14^)^ . 

^ ^j/GVo zGV^ / 

Analogously to the above, we can bound the remainder term as 

7^l(Uo, Vbn) := capf,{0,y - x) + Y capA(0, y - a;) 


xevo 

y^Vo 


xGVm 

y^Vm 


c 

< — 
- L 


sc 


revo 

y^VQ 


x-y\ L-\x-y\ 


) < C{2i+1) 


log A 


This leads to 


E E capA0. = -,)(/( " + ‘'+f 


m—l y^^o 


< = CclM log L 0 


( 6 . 21 ) 


( 6 . 22 ) 


(6.23) 


with (6.14 1 , since M < M. To conclude the proof it remains to identify the second term in the first 
line with f{u) in the limit L ^ 00 when a;/A —u G T. To this end, we use the representation 
caPA(0, z-y) = I {jY] + L-\l-v\ ) and note that 


V = V 

I r — 7/1 


2f + 1 


U-2/1 ^ ^ ^ m + f/(2f+l) 


(i+o( 


TO + f/(2£+l) 


)) . (6.24) 


Analogously, we have 
1 


E 

ye Vq 


2 £+ 1 


L — \z — y\ M — m — £/{2£+l) 


(i + o( 




M -m-£/{2£+l)JJ ■ 


(6.25) 


So the contributions to (|6.23[) of the correction terms in (|6.24|i and (|6.25|l vanish, since we have 


M 


M 


2f+l ^ 2£+l 

A {m + £/{2£+l)Y - ^ ’ 


(6.26) 
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and an analogous bound holds for the second correction term. Therefore we can replace the second term 
in the first line of ( 6.23|l by 


M 2^+l 


M 

S(; 


1 


1 


Zch L (21+1) M — m — i/(21+1) 

Rewriting the rates in this expression as 

1 1 1 

+ 


)(/( 


x+^+(2^+l)m 


L 




m+f/(2f+l) M-m- £/{2£+l) M (m/M + £/L) (1 - m/M - f/L) 


and using that f/L —^ 0 and M/M 1, (6.27 1 converges to 


(6.27) 


zJb Jo - v) 


{f(u + v)- f(u)) dv = C^f(u) 


(6.28) 


as L —>■ 00 a nd xj L —)■ w € T. By regularity of /, u i— f(u) is a uni form ly continuous function on 
T, and with (6.23 i this implies ( 6 . 81 and finishes the proof of Proposition [3^ □ 


7 Capacity estimates 


Upper and lower bounds for the capacities appearing in the proof of Proposition 3.4 in the previous 
section are obtained following closely Sections 4 and 5 in ||2l. The extension of these methods in lfT9l 
lead to matching upper and lower bounds for complementary sets of wells also in the limit L ^ oo with 
diverging particle density N/L —>■ oo. Using precise estimates on the stationary measure summarized 
in Section]^ we are able to improve estimates for the upper bound error that allow us to further extend 
these results to the the case N/L p> pc- 


To simplify notation, let N be the typical number of particles at the condensate, 

N:=N-p,L. (7.1) 

Recall that A denotes the lattice A = 'LjL'L, and we will use the subindex /c in A^ = Z/fcZ whenever 
we refer to a different lattice. 


7.1 Lower bound 

Proposition 7.1 Let Abe a non-empty subset of the lattice A and denote by the corresponding set 
of wells. Then, there exists a positive constant C such that 

iV''+^cap(£:^,£\£'^) > ^ capA(a;,y)(l-C'(L/3i"'’ + L"^+Le 2/3iV^/E)^ 

y^A 

provided that L is large enough. Here Ib = fg u^(l — u)^ du, and for any x, y £ A, capA(a;, y) denotes 
the capacity of the simple symmetric random walk on A. 

Proof. The variational formulation of the capacity establishes that cap(£’"^,£ \ £^) = 

infpg^( 5 A) 77(7^), where 

n(£^) = {F : F(p) = 1, ry e U(p) =0,p££\£^}, 
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and V is the Dirichlet form 


= ^ E E E ^^iv)9{v.)\[F{rl^’^^l-F{rl)\'■ (7-3) 

r)eXL,N xeA r= —1,1 


Define the tube 

'^L,N = Xl,n ■ Vx + Vv > N - ri:, < ISl, z ^ x,t/| (7.4) 


connecting the two wells at x and y. For any function F : > K in we can bound the 

Dirichlet form by 


V{F) > Y. 

x^A, y^A 


/ 


1 

2 


V 


E 

zGA 


riGTl’l 


(7.5) 


This holds because with Pl ^ ul the tubes only overlap within wells on which the function F is 
constant, equal to 0 or 1: DT^ C U £'^ and necessarily x = x' or y = y' for N, L large 

enough. 


Fix X £ A, y ^ A. Let Sz denote the configuration with one and only one particle at z, and for 
rj G Xj^ jy let ^ = 77 — 6z, where summation of configurations is performed componentwise. Given 
M, iT e N and C G Xm,k we will denote by p!(C) = HugAm 9KCu) = HnGAM 
rates in ( |2.3| l. Define the set of configurations = {■C G -^l, 7 V-i : ICa; + Cy ~ -^1 ^ <al/2, < 

Pl — ^ X, y} and note that ^ G ^ G for all z G A. Then we can rewrite the 

parenthesis in the right-hand side of (173]| as 


zGA SeXi,_jv_i y ts; 


r=-l, 1 f + S^GT^''" 


> 


E 

cc T^’y 


^zl.n 9m m 2 


^ E l{F{c+s.+r)-m+s.))\ 


z^A 


(7.6) 


with 


Zl,n 


E 

V^^L,N 


1 

9m 


zm{SL = N]. 


(7.7) 


If F(^ + (5a;) pG F(^ + <5y), then in full analogy with ||2l we may define 

/ = /5 : A ^ K, /(z) = 


Note that 

fm = 1, f{y) =0, and f{z + r)- f{z) = 


F{^ + ^x) — F{^ + 5y) 

F{^ + Sz+r) — F{^ + Sz) 


F{^ + Sx) — F{^ + Sy) 

so we can estimate the last sum in ( |7.6| from below by 2Lcap^(a;, y) to get 

77x.y(7^) > ^^capA(a:,y) Y {F{m ^x) - F{^ + dy)Y . 

Zl,n g'm 

L,N 


(7.8) 
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Clearly, this bound holds trivially also for terms where F{^ + Sx) = F{^ + 6y). Let be a config¬ 
uration in A \ {x, y}. For such a C G XL- 2 ,k we define the function Gq : {0,1,..., iV — fc} —>■ K by 
Gi^{i) = F{() where ( G coincides with ( outside {a;, y}, and Cx = i, Cy = ~ ^ ~ Let us 

also define the set Oi = {/c G N : \k — pcL\ < ^ ]. With this notation ([TTsIi can be written as 


'Dx,y{F) > 


Lc&-pj^{x,y) 

Zl,n 


E E 

feGOi 


1 (Gc(» + l)-gc(»))' 

ffKC) ^ 


(7.9) 


Note that if \k — pcL\ < ^ and ^ G XL- 2 ,k with max^ Cz < then = 0 when i < I3 l for L 
sufficiently large, because ^ G S'^. Likewise, GXi) = 1 when i > N — k — /3 l for sufficiently large L, 


because C, & . Hence, the rightmost sum in (7.9 1 reduces to 


i=l3L 


{N -k-i-l)^ 


Minimizing the preceding expression under the constraint 


N-k-PL-l 

E + 1) “ = 1 : 

i=PL 


the estimate in ( |7.9| l becomes 

LcapA(a;,y) 


T^x,yiF) > 


E E 


1 


k&OL<^^L-2.k 


' N-k-l-PL 

i^{N-k-l-if 

i=fiL 


(7.10) 


Using the fact that for any function / : (0,1) —>■ K, m G N and sequence km m, the following 
bound holds 


1 

m 


n—k-m 


fcr^-1/2 


i—kn 

we get that 

N-k-l-Pi, 


m/ I fcm-l/2 


fix) dx 


< 


1 


24m^ 


sup |/"(x)|, 

xGfO.l) 


Y {N-k-l-i)^ <iN-k-l)'^^+^h(\ + 


C 


i=dL 


(iV-/c-l)2 


where ly = x^(l — x)^dx. Hence, the estimate in (7.10i further becomes (provided N,L are suffi¬ 
ciently large) 




1 


> 


> 


ZLNib ^ ^ o!(C)(A^-fc-l)2f>+i 

“ keOL(^^L-2.k ^ '' 

z^-‘^Lca.p^ix,y) 

(N-k) 


( 


X 1- 


G 


Zl,n h 


(iV-fc- 1)2 

E/ (N — h\‘^b+l ^ ['^.£'-2 “ Ml-2 < Pl — f\ 


keOi 


zf ^LcapA(x,y) 
Zl,n h 


'({N - 


(7.11) 
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where Kl is the event Kl^ { \Sl-2 - PcL\ < < Pl- 1}, Ml-2 = max^gA^.a Px- Just 

as in (9.6 1 , if we define/c* = v(^Sl- 2 ', then 


+ 1^(1. -fti). 


_/\r2f)+l L J ]\J-2b+2 

In a fashion similar to ( |9.7| ) and ( |9.8| l we have 

^[^l] = v\^Ml-2 > /3l — l] + i^[|5'l-2 — PcL\ > ul, Ml-2 < h] < C (J^/3l ^ + ® 

and 

0 < PciL -2)-h<C L^PI-K 


In view of the preceding estimates (7.11) becomes 


N'^+^V^^y{F) > + +£e 


and the statement of the proposition follows from Proposition 9.1 since F was an arbitrary function 
F \ —y M, with = 1, F\^\^^a = 0. EZI 


7.2 Upper bound 

Proposition 7.2 Let A a A be a non-empty subset of the lattice, and denote by the corresponding 
set of wells. Then for any cl satisfying 


eL»^ 


and 


_ 3 ^ \A\L^N^-^ 

^ - a^^-\L-\A\) 


(7.12) 


we have 


N'‘+^ca.p{S^,8\S^) <(l + eL)j^ F, capA(a:,y) . 


(7.13) 


x^A, y^A 


In ( |7.13| ) /{, = fg u^(l — vf) du, and for any x, y G A, capA(a::, y) denotes the capacity {6.6) of the 
simple symmetric random walk on A. 


Construction of the test function. For 1 < y < L, let /o,y : A —[0,1] be the function that realises 
the capacity capA(0,?/) between the sites 0 and y for the symmetric, nearest neighbour random walk 
in A. By elementary results of potential theory, it is known that fo^y{z) equals the probability that the 
random walk reaches 0 before visiting y, when started at 0 . The formula for /o,y can be easily derived, 

r 1 - ^ if 0 < z < y; 

fo.yiz) = I " (7.14) 

[iE^y if y<z<L. 

Given e > 0, consider the smooth function : [0,1] [0,1] given by 

1 find) 

H^{t) := — / u^{l — u)^du, (7.15) 

h Jo 
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where Ib = Jg u^(l — u)^ du, and : K —)■ [0,1] is a smooth, non-decreasing function such that 
(t>e{t) + -t) = 1V< e M, <?i'e|(-oo, 3 e] = 0, <?!>eI[i- 3 e,oo) = 1, sup„gR |(m)| < 1 -f Ce and 

sup„g[g 1 ] \4't{u) — u| < C'e, for some universal constants C, C > 0. It can be easily checked that 

= i7,|(_oo.3d =0, i?.|[i-3..oo) = 1. (7.16) 

Enumerate the sites in A as 0 = Si, X 2 :..., xl = y, &o that fo,yixj) > fo^y{xj+i) for all j. 


Given rj € Xi^ jss, let fj € be given by 


fiy = Vz- Pc, 2/ e A. 


(7.17) 


For y ^ 0, define Fj : Xl,n —>■ [0,1], 1 < j < A — 1 


FlJy) = 77.(1 


N 


Po,yM " ^ I ™ ^ 


1 

N 


(7.18) 


Finally, for y 7 ^ 0, let _Fo,y • ^l,n —^ M be the convex combination 


L-l 


7^o.y(t7) = [hyi^j) - fo,yixj+i)] Fiyiv) ■ 
i=i 


(7.19) 


To define Fo,o, consider a C^, non-decreasing function /i : K —>■ [0,1] such that /i|(_oo 2 e] = 0, 
/i|(i_ 2 e,+oo] = sup,j, h'{x) < 2, and let 

Fo,o{v)--=h{yo/N). (7.20) 

In order to construct the candidate test functions for the capacities in ( |6.2| ) we need to stitch the 
functions {Fg^ylygA together. 


Let us hence define the sets S' = {it G R'^^, ^ 7 ^ 2 / £ TV, 

= {u G S : Mo + tty > 1 — e} , 

,4. = {m G S : Mo > 1 — e} 


and the disjoint sets 

jC0,y = ( j - o,y \ n {m G S, M, < |, z ^ 0, 2/} . 


(7.21) 


Let now 

{0°’" = 0°’"(i,e)}ygA (7.22) 

be a smooth partition of unity of S, such that 0°’*^ : S —>■ [0,1], t/ G A, J^y 0°’*^ — 1’ 

0°’^| o,s = 1 and sup |0 °’*^(m) — 0°’^(m')| < — ||m — m'|| 2 , (7.23) 

u,«'gs e 
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C a positive constant independent of e, L and y G A. To construct this partition of unity, consider 
a bump function p : M — >• [0, 1], p = 1 on (—oo, |e], p = 0 on [e,+oo), |ip'||oo < Let 
d{u, = inf{||u — w\\ 2 ,w G K be the Euclidean distance from a point u to the set 

and let 0^’^(u) = g(d(u, IC^’^)), 1 < y < L — 1. Note that these have disjoint supports and take 
values in [0,1]. Dehne = 1 — X)i<j/<l-i 

The candidate to solve the variational problem for the capacity cap(£’°, £’\£°) is then Fq : at —>■ 

M, 


f'ojg)- (7.24) 

yeA 

Lemma 7.3 The function Fq in ( |7.24| ) satisfies 

Foir,) = Fojn), forfjGlC°’y, (7.25) 

Fair,) = 1, for % > (1 - 2e)N, (7.26) 

Fo(77)=0, forfio<2eN, (7.27) 

\Foiv)-Fo(v')\<C^ \\ri-T]'\\ 2 , far rj, g'€ Xl,n. (7.28) 


Proof. The hrst assertion follows from ( |7.23| l. For the second one, note that = 1 on [1 — 3e, 1] hence 
Fo,y(v) = Fo^y(r]) = 1 for all ?/ G A \ {0}, j = 1... L - 1, and also Fafag) = 1, if ijo > (1 - 2e)N. 
Similarly, ( |7.27| l follows from F^^fa) = Fo^y{r]) = 0 for all y G A\ {0}, j = 1... L — 1, To.o(??) = 0, 
if Vo < 2eN. Finally, ( |7.28| l is a consequence of ( |7.23| ), if we use that the supports of the functions 
00-^ y f 0 are disjoint, and hence at most four terms in the difference 

Foil) - Fq[v') = E FQ,y[v) - 0°'\v'/N) Fo,,(p')] 

yeA 


do not vanish identically. □ 

By \1 .26) and ([1 .21) , Fb = 1, Fq hence we can use it to estimate the capacity between 

£0 and £ \ £°. 

Single well capacities. Given / : at — > M and U C dehne 

CNifaU):=lY.Y. E T{v)9ivfal[fiv^’^+l-f(v)]\ (7-29) 

r}GU xGA r—— 1,+1 

and for x, y G A a let 

■= {v G Xl,n, Vx+Vy>N- aL}, (7.30) 

where s/~L az, A is the constant used in the dehnition of the wells ( |2.12| ). 

We hrst notice that the only relevant contributions to the Dirichlet form ( |7.3| l originate from conhg- 
urations in the sets y G A. 
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Lemma 7.4 


Cn{Fo; \ Uo^yeA^jv^) — ^ ' 
C a constant independent of L and e. 




2 AT2 2b—2 ’ 

r ar 


Proof. By ^TXT) , ( IXlSj l in Lemma[73 


< C 


L 

g2JV2 

L L 


Cn{Fo',Xl,n \ Uo^ygAlJ^^) < C P & Xl,n \ Uo^yeA^^^; 

L2 


2e < 3 < 1 - 2e 
~ N ~ 


}] 


= C 


2 AT2 2b—2 


The last line follows from Proposition |9.6| in Section]^ 
Note that 

Civ(Fo;Uo^y6Ai^^) < 

yGA 

as all terms in the latter sum are non-negative. 


□ 


(7.31) 


Lemma 7.5 Let e > 0 be such that eN ^ a^- There exist positive constants C, C which are indepen¬ 
dent of L, and such that for any 0 f y G A, 


\Cn{Fo;I°^^) - C^(Fo,y;I^*')| < 


(7.32) 


Proof. By Lemma 7.3 we may just consider those configurations p such that 2eN < Pq < {1 — 2e)N. 
Note that the fact that p € implies that fjy > eN, and also that for Fo{p) f- Fo^y{p) it is necessary 
that 0°’^(7)/iV) f 1, which combined with the previous cond ition s implies the existence of z G A, z 
0, y, such that pz > eN /2. From ( |7.29| l and ( |7.28| l in Lemma 

|C;v(Fo;I^^)-C'^(Fo,y;Z^*')| 


7.3 


we get 


<C-^p 

- £2^2 ^ 

<C^^p 

- e2jY2^ 


{n e 2e < ^ < (1 - 2e); % > e; 3z ^ 0, y, i);, > ^ 


N 

L-2 


{3=^0.!,, y: 


i'>2,Xi^z 


(7.33) 


By Proposition |9^ 


\Cn{Fo-,X°m^) - CN{Fo,y-,I^^^) \<C-^ exp{-C'e2iVV(L - 3)}. 


with 




> 


L-3 — L ■ 


□ 


Proposition 7.6 Consider e > 0 such that eN ^ aL- Let y f 0 in A. There exists a constant C > 0 
independent of L, e and y, such that 


C^(Fo,y,Z^^) < (1 + Ce) 


Lz, 


L-2 


capA(0,t/) 


7V2''+1 Zl,n h 

capA(a;, y) the capacity between sites x and y for the simple symmetric random walk in A. 


(7.34) 
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Proof. We first focus on terms r = 1 in \1 .29) . Fix i S A, and let 1 < *, j < L be the indexes of i, i+1 
in the enumeration of A determined by sites 0, y. Let us suppose i < j, that is, fo,y{x + 1) < fo,y{x). 
By ( |7.19| l, we get 


f-i 


- FM = E - hvi^k+i)] [Fly{ri^’^+^) - F^o"v(7)] • (7-35) 


k—i 


Replacing ( |7.35| l in definition \1 .29) , 

E ^^i^)9{vs:) I [Fo.yiv'"'''^^) - Fo.yir])]' 


veil 


f-i 


— [fo,y{x) fo,y{x + 1)] \fo,y{xk) fo,y{xk+l)\ 


k—i 


veili'> 


(7.36) 


by Jensen’s inequality. Note that terms associated to p € and po < 2eiV or 770 > (1 — 2e)N vanish 
by the definition of F^y, ( |7.18| l. 

Let 

i 

:= {77 e 2eN < 770 < (1 - 2e)iV, 3 2 < ^ < L - 1, | ^ 77 ,, | > eiv}. (7.37) 

i^2 

Fix now i < k < j — 1 and consider the sum in the last line of ( |7.36| t associated to the index k. We have 


7-0)y\ >aO>y 
\^N 


+ E^( 7 )^ - PoJv)]' 


where 


g0,» ^ 

^ l\T 


(7.38) 


In we have < a^. Together with the fact that there is 2 < Z < L — 1 such that 

I J2i=2 V=^i\ ^ implies that either a) Yl!i =2 ^ Yld=iXi Vxi < ocl - eN- Once 

more, by Proposition [E] and with the hypothesis that eN ^ a, we get 

< Cl(7.39) 

Cl, C 2 positive constants that do not depend on L. From ( |7.39| l and ( |7.38[ l it follows that Civ (Ab,y > 
is negligible to any polynomial order in (|7.34|l. 
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On the other hand. 




7-u,y\ wO 
-^iv \^i\ 


< 


E 




-. K ^ AC ^ 

rl+E^-i/.fEM 

\N '^TV/ V^iV/ 

r=0 "— 


with g!(C) = no<i<L-i 5’(6) = no<i<L-i where J°ti is the set 


r—0 


(7.40) 


Jn -1 ■= {ee^L,^-i;2e< I < il-2e), io + iv > N-aL;\^U\ < eNy2 <l < L-l}. 


By the dehnition ( |7.15| ) of H^, for ^ G we have 


1 ^ t7 ^ <7 

HJl + y^kA-uJy^kA 

\N ^ TV / N / 


^ TV 

r=l 


r—1 


^+Eo<r<fc 


4^ 


.E 


(1 — u)^ du , 


0<r<fc jv y 


and hence 


. „ 0 ,» i' ''S’/' L r=0 r=0 




< 


gm 

Qi —1 

E 


E 


m——pc{L—2) ^^Xl,-2,m ■, 

M—pc{L—2)-\-m 


1 1 1 


(1 —2e)JV+Pc / 

^ i: U 

Z.=2€iV + Pc ^ 


u^(l — u)^du 




(7.41) 


where 4>, := 


:= + i lair)- order to derive the right hand side in ( |7.41| l, we consider 

the distribution of m + Pc{L — 2) particles in the L — 2 sites other than 0 and y, and then the 
distribution of the remaining particles between these two sites, such that t particles are assigned to 0 
and TV — TO + 2pc — l — 1 particles to y. 


Terms in ( |7.41| l associated to indices —pc{L — 2) < to < — eTV correspond to moderate or large 
deviations events 'Yli 2 <i<L ^xi — Proposition 9.7 implies that their combined contribution decays 
faster than any power of L, and is hence negligible. 

Now, if —eTV < m < and recalling that | X]i=2 4;; I E eN^l, we have 


-r- ^ du < ($ 1+1 - ^ 

(t)fa(l m+t-2p^ + l \b — V i+i 


l-$, 


N 


N 


-| m-j-A-—2pc + l 

. N / 


< ($ 1+1 - $t) sup 

u£ [2e,l —2€] 

< i (l + C'e)2^ 

~ N 


Mu) 

u- fe 


2b 
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In order to obtain the last line, we applied Cauchy’s mean value theorem to ^ using 


that ^e(3e/2) = 0 and sup„ \(j)'^psilon{u)\ < 1 + Ce from the dehnition after (7.15 i. Replacing this 
bound in (|7.41|l yields 


. ^o,y |_ r=0 r—0 




<(1 + Ce) 5] 


E 


1 1 1 (l-2e)tV+p„ 

E / 


1=-Pa{L-2) ieXL-2,M, 6 _ ^_2£Af + Pc 


N^b+1 II g\{0 


u^(l — u)^ du 


M—pc{L—2)-\-m 


i=2 


< (1 + Ce) 


E E 


r ^26+1 o!(£) 

^ m=-p,(L-2) 4 GXz,_2.m, ^ 

M—pc{L — 2)-\-m 


^ < (1 + Ce) 


L-2 


4 Af26+1 


We can now put the pieces together: replace the above estimate in (|7.40|) to derive from (|7.36[) that 


E^(,) < (1 + c.) — 

'n^^N 


[/o.y(®)-/o,y(5 + l)]' 


We conclude that 


CN{Fo,y,I°^^) < (1 + Ce) 


^L-2 


= (1 + Ce) 


r^E E ^[/o.y{*)-/o.!/(S + k‘ 


Af2b+i Zln Ib 2 “ — 

x£A r= —1, + 1 

capA(0,?/) 


iV2b+i Zl,n h 

The extra factor L in the line above is compensating for the uniform weight ^ appearing in the 


computation of the capacities for the symmetric simple random walk on the torus. 


□ 


Lemmas |7.4| |7.5| Propositions |7.6[ [9T| and observation ( |7.31| l combined yield an upper bound for 
the capacity 

N’’+^cap{£°,£\£°) < -—^capA(0,y)(l + (Ce + - + --). (7.42) 

V L e^a^ logL/ 

By the symmetry of the model, a similar estimate holds for cap{£^,£ \ £^), x G A, evaluating in 
this case the Dirichlet form at = Fq{t-x(,v))’ if '''-a; ■ A^l,n —>■ Xl,n is the shift [T_a:(? 7)]2 := 
z G A. 


Multiple wells capacitites. As described in Section]^ to compute the rescaled rates of the process we 
need bounds on the capacities 

cap(£:^,£\£:^), (7.43) 

where we recall the notation £^ = L)z(=a£^, A C A. The strategy is to compute these in terms of the 
single well capacities estimated in the previous paragraph. 
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At this point one would like to take advantage of the subadditivity of the capacities to bound 
cap(£'^,£ \ £^) < ExeAMA cap(£®,£^). To make it work, instead of estimating the capacity 
between complementary sets cap(f“, £ \ £^) in terms of capA(a;, y), as we did in the previous 

paragraph, it would be necessary to estimate each term cap(f“,£^) directly from the associated 
capacity capA(a;, y) for the simple symmetric random walk on A. Unfortunately we are not able to do 
this: we need the function to be the relevant term on the set but it also contributes to the 
other sets , z ^ y, in non-negligible ways. Indeed, it seems that the only solution to this problem 
is to consider a partition of unity that singles out each of the sets z ^ x, as in ( |7.21| l, and to 
dehne as in ( |7.24| i, thus reducing the Dirichlet form to the sum of the Dirichlet form restricted to 
these sets, while still providing a suitable test function 77^ such that F^\g^ = 1, Fx\£\^g^ = 0; that is, 
the construction of the previous paragraph. Following Q, we next show how these functions can be 
combined to estimate cap(£"^, £ \ £^). 


Proof of Proposition 

EzgA SO that F^ 


7.2[ Fix Cl > 0 be as in the statement of the proposition. Dehne F^ := 
= 1 and F^ L,^ = 0 by ( |7.26| l and ( |7.27| ). As a hrst observation, no- 




tice that it is enough to consider 

Indeed, using Cauchy Schwarz, 

\ y^Z ^ ^ ^ ^ 


(7.44) 


< C- 


zGA 

\A\^L^ 




El 




by Lemma 


(7.45) 


ce Xl,n \ Uzga ( UyeA C Xl,n \ Fiy^u for each given u G A. 

Next, 


Cn(f^; U,eA( U,eAX^^ )) < Cat 

= E + EE+ E 


z^Ay^A 


u)/z, y 
wGA 


zeA y^A 


w^z 

wGA 


(7.46) 


We will consider the terms in these sums individually. 


We have 


w^z, y 
w^A 


CN^F^+FyF Y, F^;r^yyCN[F,+Fy;r^y^\<C^^ti[F^^yn{fix>2eLN,x^y,z} 




on account of ( |7.27| i and ( |7.28| l. Now on n {fj^ > 2eAiV} we have J2u^z y x^u — ^ ~ ~ 

oil) — ^lN = aL — 21 lN, and by Proposition |9^ 

n {nx > FIlN, X + y,z}] < 


Cl, C 2 positive constants that do not depend on L. Furthermore, the arguments in Lemma 7.5 also yield 

|Civ(x, + Fy; X^*') - C^y^y + Fyy, X^«) | < Cie-^=“"/^, 
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so that 


CN{F, + Fy+ FfX/) 

w^z, y 


CM{F.,y+Fy^,-r^y 




z,y&A, (7.47) 


where the values of the constants may change from line to line. By similar estimates, 

\Cn[f,+ ^ z G A,y ^ A. (7.48) 

w^z,y 


We next show that terms C^i^Fz^y + Fy^z, F^^j, y,z G A, are negligible to first order. Let 77 G 
\ U (see definition ( |7.37| i). We get 


FUy) + Fy^r^y) = ff,. 





= H-,. 







= 1 . 


The first equality holds by the fact that the enumerations and {xj}^ ^ on A determined by 

fz^y and fy z satisfy Xi = The second equality is due to the identity (f) + (1 — f) = 1. 

Then 


Fz,y{y) + Fy^z{y) = 

L-1 L-1 

= [fzA^j) - fz,y{Xj+i)] Fiyij]) + ^ [fy,z{Xj) “ fy,z{^j+l)] 

i=i i=i 

L-1 L-1 

= Y “ L.y(^i+i)] + Y - h,vi^j+i)] 

i=i i=i 

L-1 

= Y [L.y(a;j) - fz,y{^j+i)] {F'lyiv) + F^-^{y)^ 

i=i 

L-1 

= Y - h,vi^ 3 +i)] = 1 > yGl^jf\ {B^jf U B^/), 

i=i 

where we used that fy,z{xL-j) - fy,z{xL-j+i) = fz,yixj) - fz,y(xj+i). In particular 


Cn {Fz,y + Fy^z , \ U = 0. (7.49) 

Finally, the arguments leading to ( |7.39| ) yield 

Cn [Fz,y + Fy^z , B]f U _ (7.5O) 

Collecting all results from equations ( |7.45| l to ( |7.50| l, we obtain 

Cn{.F^-. Xl,n) <YY + 

zGAy^A 
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and by Proposition |T^ 


Cn(f^,Xl,n) < (1 + Ccl) j^2b+lzr „ 4 ^ ^ capA(0,2/) +t^-2^72„>-2 ' 


L-2 


|A|2L2 


^L,N 




In particular, pulling out the first term as a common factor, and applying the rough estimates 
capA(z, y)> ^ and X;^gA Ey^A capA( 2 , y) > condition \1.12) yields 


iV''+icap(£’"^,£:\£’"^) < (l + Cci) 


Lz, 


L-2 


Zl M Ib 


^ capA(z,2/). 


zGA, y^A 


The assertion of the proposition follows from Proposition|9.l[ as eLL ^ 1/L. 


□ 


8 Tightness and Limiting Distribution 

8.1 Proof of Proposition Tightness 

In this section we will prove tightness of the distributions of (Y/" : t>0) := : 

f > O). By Aldous’s tightness criterion (cf. Theorem 16.10 in ifTTl l. it suffices to show that for any 
77 S £, e > 0, f > 0, 

limlimsupsup sup Fr,[dj(Y^^g,Y^) > e] =0, (8.1) 

i54^0 L —>.00 s<S 

with dj{x,y) = \x — y\{\ — \x — y\) the distance in the torus T, where, as before, all algebraic 
operations are performed modulo Z, i.e. x — y — (x — 2/)mod(z)- ( |8.1| ) Tj is the set of stopping 

times (for the trace process) bounded by t. 

Note that by the strong Markov property, for any y G S we have 

P, > e] = E, [P,.(.) [driY,^, Y,^) > e]] 

< sup^g£ Py [dT{YL Yq^) > e] • (8-2) 

We will denote by 

r^{v,z) 

xeA 


the total jump rate from a configuration 77 in some well to a new well By Ito’s formula, for 

any s > 0 we have 


dj{YLYo^) = 


rsOr 


{rj^(u); z) di(V'i(77(77)) + z/L, Fo^) - dT{'iljL{v{u)),YQ^) 


z#0 


du + MsBl 


where {Ms}s>o is a martingale with quadratic variation 


{M)s = / ^ ( 77 ( 77 ); z) [dT(7/'L(7?^(7i)) + -dT(7/’L(77(u)),Fo^ 


du. 


47 













By the triangle inequality we get 

IdTitpLiv) + zlL,Y^) - dT{'iljL{v),yo^)\ < dT{z/L,0) = 
By Doob’s inequality, for any rj G £, 


■ V 


1 C p 1 1 C _ 7^ / 7\2 r 

sup \Ms\>- <^E, {M)se, <^^.-^[^- 7 ) / 

^sKser. 2J L J e ^ LJ IJq 


16 




z^O 


pse, 


[t]^{u); z) du 


(m); z^du 

(8.3) 

(8.4) 


On the other hand, by Markov’s inequality, for any rj G £, 

P,,[ sup f y)r^{r]^{u);z)du> ^ 

'-s<S0r. Jo ,r, ^ Z 


z^O 


“ S<S0L Jo 

2^0 

Hence, it suffices to show that 

lim lim sup sup — (1 — — ) E. 
<510 r_>oo r,GE,^L\ Ll 


n56 T 


‘-^0 


[t]^ (u)] z')du 


2#0 


n59 T 




[■q^{u)] z) du 


= 0 . 


(8.5) 


( 8 . 6 ) 


We proceed analogously to the proof of Proposition |3.3| giv en in Section |4T| using the bounds on the 
mixing time fniix(eO for restricted process given in ( |4.3| l. We pick e' = 1/6^ and split the integral 
in the preceding display according to 


pS0, 


(rj^ (u)] z) du = I (r]^ {u); z) du + I {rj^ (u); z) du . 




p69r 


^mix ) 


We can estimate the first contribution to ( |8.61 l as 

/*^mix(^ ) 


25^0 


1 - - E 


^Jo 


(r]^(u)]z) du 


f mix(e') 

4 


< sup ^ / (p; z ). 

2/0 


(8.7) 


( 8 . 8 ) 


By translation invariance sup^g£ J2z^o = sup,jg£o J2z^o (7> which vanishes as L ^ 

00 by Lemm as 1 3 .5 1 and |4~T| For the second contribution in ( |8.6| l we use the definition of the mixing time 
and the fact ( |4.21 l that the invariant measure fY of the restricted process equals the invariant measure fi 
restricted to the well f to get the upper bound 

n59 T 


EK 


4 1 -^ E 


2#0 


^mix ) 


(r]^{u),z) du 


^ SdLe' 

- 

4 neso 


Y^r^(rj,£^)+6e,Y.y[l-y)r\z). 

2^0 2#0 


The hrst term vanishes as L —00 with the choice of e' above and Lemma 3^ and for the second term 
( |6.8| l implies 


2/0 


zjbip - 


hence establishing that the family {Q^}lgn is tight. 


(8.9) 

□ 
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8.2 Proof of Proposition |3^ - Martingale convergence 


The rescaled position of the condensate {Y^ : f > O) is a random variable that takes values on the 
space -D([0, T]; T) of cMlag paths on the torus T, and as we proved in section 
corresponding di strib utions {Q^}lgn is tight. In this section we prove Lemma 
to show that in fact 


we use in section 


with a preliminary result. 


3.1 


8.1 


8.2 


the family of the 
stated below, that 


is convergent and to characterise its limit. We begin 


Lemma 8.1 Let Q be any subsequential weak limit of Q^, U a bounded continuous function in 
Zl([0, T]; T) and f : T —>■ K a continuous function on the toms. Then, for any t > 0 


U{uj)f{ujt) dQ^{uj) 


U{uj)f{u}t) dQ{uj). 


Proof. Consider the mapping : D(^[0,T]; T) —)• M with = /(w*). Then 11^’* is continuous 

at paths oj G D(\0,T]; T) that are continuous at t, i.e. uit = oJt- '■= lims^tWs. To see this, note that 
for any u G D(\0, T]; T) and At G [0, T] we have 

|n/'‘(u) - n/'‘(a;)| < |/(ut) - fM\ + \fM - /(a.t)|. 

Since / is continuous on the torus and uj is continuous at t, the right hand side of the preceding display 
can be made arbitrarily small, provided we can control | rtf — \ and | Af — f |. If we choose u sufficiently 

close to UJ in the Skorokhod topology we may find a At that simultaneously makes these quantities 
suitably small. Since the set of discontinuities of the function uj i —U (yj) f {ut) is contained in the set 
of paths {w : uJt f Wt- }, the assertion of the Lemma will follow if we show that 

Q[a;f ^a;f_] =0. (8.10) 

For f > 0, and £, (5 > 0 define the subset C f2([0, T]; T) by 

j/’® = {w e L)([0,T]; T) : sup djiujt^ujs) > e}- 

We have {ujt uJt-} = Ue>o n 5 >o Hence 

Q[a;f f ujt-] = 0 limQ[j/’®] = 0, Ve > 0. 


Noting that are open in the Skorokhod topology it would suffice to show that 


limlimsupQ^ [j/’®] = limlimsup P[ sup djiY/",Y^) > e] = 0, Vi5 > 0. 

i5i0 too L—tcjo 


But this follows from tightness estimate in Section 8.1 in particular (8.3 1 and (8.51 and estimates after 
that. □ 


Note that for any Lipschitz function / : T —> K we have that £^(/) is a continuous function on T. 
Indeed, using the elementary estimate 

\f{y + u) - f{y) - f{x + u) + f{x)\ < 2Up{f){dT{x,y) A dT(M,0)) , 

we have that 

\^'^f{y)-d:'^f{x)\ < 2Lip(/)^ ^1 A du < AY\p{f)dr{x,y) ^2 + In . 
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In particular, the mapping w i— C'^f{uJs)ds is continuous in the Skorokhod topology, and combined 
with Lemma 8.1 we get the following; 


Lemma 8.2 Let Q be any subsequential weak limit of and / : T —>■ M a Lipschitz continuous 
function on the torus. Set 

Ml (w) = f{ujt) - f{uJo) - [ C^fiuJs) ds. 

Jo 

IfU is a bounded continuous function in I?([0, T]; T) and t > 0 we have 


j U{uj)mI (w) dQ^{uj) 


U{u})mI (w) dQ{(jj). 


8.3 Uniqueness for the martingale problem 

In this subsection we prove uniqueness for the martingale problem associated with the operator 
i.e. there exists a unique measure Q on the space Z3(IR_|_;T) of cMlag paths on T, such that for the 
coordinate process (ujs : s > 0) we have 


= 2/0 e T] = 1, 

f{oJt) — /g jO.'^f{oJs) ds is a Q — martingale for all / £ Lip(T). 


( 8 . 11 ) 


We may plug the test function fk{x) = x £ [0,1], fc € Z in (8.11 1 to get that if Q solves (8.11 1 , 

then for any f > s > 0 we have 


I Ts] = where -f{k) = H{b, p) [ 

Jt 


1 — cos(27rfc?/) 
c(t(0,2/) 


dy. 


In particular, this shows that under Q the coordinate process has independent, time homogeneous incre¬ 
ments. This determines the hnite dimensional distributions of (ojs '■ s > 0) and implies that such a Q is 
unique. To get a better insight on the limiting process, we may rewrite ip as 


iPik)=H{b,p) 


1 — cos(27rfcy) 

\y\ 


dy 


and note that if (Xj : f > O) is a symmetric Levy process on K with Levy measure 

\y\ 

then Q is the distribution of the process (Xt(mod 1) : t >0). 


9 Estimates on the invariant measure 

In this section we collect some auxiliary results of technical nature that we needed throughout this arti¬ 
cle. We will use C to denote a constant (not always the same) that only depends on absolute constants. 
When X is a random variable dehned in a probability space (12, X, v) and X € X we will write v{X] 
as a shorthand for X dv. Recall also that N = N — p^L, and is the normalization in ( |7.7| l. 
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Proposition 9.1 Suppose S > pc- Then, provided that L is sufficiently large, we have 


sup 

N>SL 


Zciy[SL = 


LN 


-b 


- 1 


c 

< — and thus sup 

L N>5L 


Zl, 


N 


z^-^LN-b 


- 1 


C 

< — . 
- L 


(9.1) 


Proof. The proof essentially follows the argument in IISl keeping track of the rate of convergence 
when N is supercritical. For a given sequence we define the events = {rjx < h^}, and Ck = 
n • • • n n Bk+i n • • • n b^. Then, 


[SL = N]=j^{^\[{SL = N}f^Ck\. 

fc=0 ^ ^ 


(9.2) 


Let us denote by G the distribution function of the one-site marginal of the critical measure, i.e. G{t) = 
b'\rix < f] and let G = 1 — G. If is chosen so that LG^Hl) 0 hi, ^ , then Lemmas 2.4 

and 6.1 in ll23l together imply that 


Y^^(^fj^[{SL=N}nGk] <CLi,[{SL=N}nGi] X (lG(/il)). 

We will choose hi = so that the sum of the terms for fc > 2 in (9.2 1 is at most 0{L~^) times 

the term for k = 1. To estimate the term for fc = 0 we may use Lemma 2.1 in ll23l . Precisely, if 
K = max{ 2}, then L” is a natural scale for Si — pffi and 

v[ \Sl- pcL\ > X, Ml <hL\ <v[\S l- pcL\ > x, Mi < L''] < Ge“^^ " Va; > 0. (9.3) 

Since N > {S — Pc)L and k < 1, for all sufficiently large L we have " < L~^, hence 

L i^[{Sl =N}n Gi] <4Sl = N]<L iy[{SL = N} n Gi] (l + ^) ■ 

We now pick a sequence ji = L'’', with k < 7 < 1. Using ( |9.3| l again, it is not hard to see that 

N-hx 

iy[{SL = W} n Gi] = - ^ (iV - fc)-' = k, Ml_i < hi] 


^ fc=0 

= - Si-i)-^; {5l_i <N-hL, Ml_i < hi}) 

= 1^((7V-5l_i)-'’; Ki) (l + o(y)), 
where Ki = { [Bl-i — PcL\ < 7l, Mi-i < hi}. It also follows that 

-v{{N - Si_i)-^- Ki)<v[Si = N] < -v{{N-Si_i)-^-, Ki) (l + y). 

Zc Zc L 

If IBl-i — pcL\ < ji we have the following pointwise inequality. 

0 < (IV - Si-i)-’’ - N-^ - - pffi) < G7V-(''+2)(Bl_i - pffif. 

Integrating over Ki and setting fc* = Ki) = [L — l)v{rii\ Ki) we get 

i^[Ki\ b{k^ - PcLv[Ki\) 


(9.4) 


(9.5) 


0<i/((iV-BL-i)-'; Ki)- 


Nb 


Nb+l 


G GL 

< ^— -iy{{Si_i - PcL)^) < 


]S[b+2 


]S[b+2 


(9.6) 
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We now have to estimate how close i'[Kl\ is to 1, and fc* to pcL. Using (9.3 i again we have 


i'[Kl] = v[Ml-i > /il] +i^[ \Sl-i - PcL\ > 7 l, Ml-i < hi] 

<{L- > hi\+v[ |5'l_i - PcL\ > ji, Mi_i < L”] 

< + (9.7) 

We also have 


Pc = '^{m) > Kl) ( = K/{L - 1 )) 

= v{rii-,Mi_i < hi) - 1/(771; - pcL\ > 7 l, Ml_i < hi) 

> G{hi)^~^v(pi-,pi < hi) - 7/(771; \Si-i - pcL\ > 7 l, Mi^i < L'^) 

> (1 - LG{hi)){p, - Ghl-^) - 7/(77^) > pe(l - GLh]-’^). 


In the penultimate step we have used Cauchy-Schwarz, (9.3 1 , and the elementary inequality (1 — x)" > 
1 — nx. Hence, 


2.1-6 


0 < Pc{L-l) -h< CL^h] 


Plugging this estimate into (|9.6|), the assertion follows from (|9.5|l, since Lp'h]^ ^ = 1. 


(9.8) 

□ 


In the proof of Proposition |9.3| we will need the following lemma, which holds in particular for 
subextensive N. 


Lemma 9.2 Suppose L ^ 00 and N ^ L log L. Then, 

Zi^ji = z]-^LN-^{l + o(l)). 

Proof. This follows immediately from Theorem 2.1 in ll23l since 

1 


^gn= n 

V^^L,N xGA 


g'-ivx) 


= z]i^[Si = N]. 


(9.9) 


□ 


Corollary 1 in 13 states that for 6 > 3 the fluctuations of the condensate size around N are of order s/L 
and asymptotically normal. The next Proposition provides a conditional large deviation upper bound for 
the size of the condensate. 

Proposition 9.3 Suppose L'^ < ai < N for some 7 > i. Then, for L sufficiently large we have 

p[Mi < N - ai] < GLa]-^. (9.10) 

Proof. We first observe that there exists a positive e = e(5, 7 ) such that for L sufficiently large we have 
p[Mi<eN]<GLN^-\ (9.11) 

This follows from Lemma 5 in ll24l and Lemma [9^ since 

p[Mi < eN] = -y^iy[Si = N, Mi < eN] < 

^L,N ^L,N eN^ 
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(L- 1 ) 0 , 


L-2 


N — OiL 


Zl 


N 


(1+ «(!)) E 


1 


(N -kf' 

k^eN ^ ^ 


(9.12) 


Choosing e < ( 26 - 1 ) 7-1 (9.11 1 holds for L large enough. With this observation the assertion 

follows easily when > (1 — e)N. Indeed, in this case 

fi[ML<N- ul] <h[Ml< eN] < CLN^-^ < CLa]-\ 

If on the other hand ul < (1 — e)A^we have 

Ijl\Ml < N — cti] < /r[Mi < elV] + Lfj,[eN < Ml = "ni < N — cil] . 

Now, 

^[eiV < Ml = m < N - ul] = — E ub < k, Sl-i = N - k] 

^L,N ^ 

k^eN 

L-1 N-oll . 

' k=eN 


(9.13) 


where the last line is a consequence of Theorem 2.1 in Il23]l and the fact that N — k > pcL + if 
k < N — aL- By Jensen’s inequality for the convex function x 1—)■ a;“^(l — x)~^, x G (0,1) we have 


N—ax, 

E 


1 


fi- 




nN 


k^{N -kf 


< N 


1-26 


dx 


T-b„.l-b 


^ < CTV-V^ , 

X^{l — xf 


for L sufficiently large. Applying this estimate to the last line of ( |9.13| l, we conclude that 

T fj-b 

p[eN < Ml =m<N- a^] < 

^L,N 


In view of Lemma 


9.2 


the assertion follows by replacing (|9.11|l and (|9.14|l in (|9.12|l. 


(9.14) 

□ 


E 


Zl,n ^ g'-(k) 

k'yN — oll 


v[Sl-i = N - k, k> Ml-1 > Pl] 


< 


Lz, 


L-l 


< 


Zl,n{N — aLf 
Lz^-^ 
Zl,n{N — ul) 


v[Sl-i < PcL + oiL, Ml- 1 > h] 
^v\Ml-i > I3l\ ■ 


(9.15) 


1 

Proposition 9.4 Suppose '< chl < for some 7 > 1/2. Let M/ ' stand for the second largest 
component 0 /( 771 ,..., 77 ^). Then for L sufficiently large we have 

p[Ml>N- aL, Mf) > h] < CLUI-K 

Proof. We have 

p[Ml > N - aL, > Pl] < Lp[pL > N - aL, Pl > Ml-i > h] 

Lzf:-^ ^ 1 
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The assertion now follows from Lemma [9^ and the elementary estimate 

y[ML-i> I3 l\<{L-1)G{Pl)<CL^I-\ □ 

Propositions |9.3| and |9.4| show that the invariant measure fjL is essentially supported in the union of 
the wells, as the following corollary states. 


Corollary 9.5 Consider the wells defined in {2.12\. Suppose < aL < for some 7 > 1/2, 
and fii, ^ . Then, for the complement A = X]^ jq \ have 

p[A] < CL{^a]f^ + 0 as L ^ 00 . 


Proof. It suffices to note that 

p[A] = p[{Ml <N-aL}U {Mf ^ > Pl}] 

= ti[ML<N- ai] +tJi[ML>N- oiL, iwf ^ > Pl] ■ 


□ 


Proposition 9.6 Suppose for some 7 > 1/2, and consider the sets defined in 

( |7.30| l. Then 


T 


{ri e \ U ■■ 2 aL<Vo<N- 20 ^} 

O^yGA 


CL 

— „2b-2 ■ 


Proof. We have 

T 


[t] e Xl,n \ U ■■ 2aL < Vo < N - 2 ai} 

O^yGA 
N—2oll ^ 

= X! ^^ ~|7TT TL-l,N-k[ML-l < N - k - aL\ 

gm 


CL 

- ^ kHN-kf 

CL dx ^ CL 

~ {aLNf-^ xf>il-xy - af-2‘ 


The third line follows from the second by Lemma [9^ and an application of Proposition |9. 3 1 term by 
term, with supercritical particle number N — k > p^L + 2aL, while the fourth line follows from the 
third by Jensen’s inequality. □ 


The following Proposition states that deviations of occupation numbers below their typical value 
have exponentially decaying probability. 


Proposition 9.7 Suppose N > L'^ for some 7 > |, and let A C X Let ijx = Vx ~ Pc v = 2 i/(t]x). 
For any m > 0 and large enough L we have 


< -m 

x^A 


< Ce~^. 
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Proof. Write 


>‘[T. 

x^A 


fjx < -m 


^ E ^[S\A\ = k]v[SL-\A\ = N-k]. 

^L.N , ^ I ,t I 

k<pc\A\—m 


For k in the range of summation the number of particles in A \ A is supercritical. Indeed, N — k — 
Pc{L — |A|) > N + m > N. Hence, v[Sl-\a\ = N — k] < C{L — |A|)7V“^ for L large enough, and 


' x£A 


lx < -m 


< 


Cz^ {L-\A\) 

Zl,n 


E ^[^\A\=k] <Cv^^ix<-m 

k<pc\A\—m xGA 


where the last inequality follows from Lemma |9.2| The rest of the proof is now a standard Chernoff 
bound for the sum of the independent variables {ix}x&A- For any A < 0 we have 


[ E - 

x^A 



= (1 + Ape + . 


1^1 


The assertion now follows by choosing A = —2m/u|A|. 


□ 
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